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ABSTRACT 

Millisecond pulsars are mainly characterised by their spin periods, B-helds and masses 
- quantities which are largely affected by previous interactions with a companion star 
in a binary system. In this paper, we investigate the formation mechanism of mil- 
lisecond pulsars by considering the pulsar recycling process in both intermediate-mass 
X-ray binaries (IMXBs) and low- mass X-ray binaries (LMXBs). The IMXBs mainly 
lead to the formation of binary millisecond pulsars with a massive carbon-oxygen 
(CO) or an oxygen-neon-magnesium white dwarf (ONeMg WD) companion, whereas 
the LMXBs form recycled pulsars with a helium white dwarf (He WD) companion. 
We discuss the accretion physics leading to the spin-up line in the PP-diagram and 
demonstrate that such a line cannot be uniquely defined. We derive a simple expression 
for the amount of accreted mass needed for any given pulsar to achieve its equilib- 
rium spin and apply this to explain the observed differences of the spin distributions 
of recycled pulsars with different types of companions. From numerical calculations 
we present further evidence for significant loss of rotational energy in accreting X-ray 
millisecond pulsars in LMXBs during the Roche-lobe decoupling phase (Tauris 2012) 
and demonstrate that the same effect is negligible in IMXBs. We examine the recy- 
cling of pulsars with CO WD companions via Case BB Roche-lobe overflow (RLO) 
of naked helium stars in post common envelope binaries. We find that such pulsars 
typically accrete of the order 0.002 — 0.007 Mq which is just about sufficient to explain 
their observed spin periods. We introduce isochrones of radio millisecond pulsars in the 
PP-diagram to follow their spin evolution and discuss their true ages from comparison 
with observations. Finally, we apply our results of the spin-up process to complement 
our investigation of the massive pulsar PSR J1614— 2230 from Paper I, confirming 
that this system formed via stable Case A RLO in an IMXB and enabling us to put 
new constraints on the birth masses of a number of recycled pulsars. 

Key words: stars: neutron - white dwarfs - stars: rotation - X-rays: binaries - pulsars: 
general - pulsars: individual: PSR J1614— 2230 



1 INTRODUCTION 

Binary millisecond pulsars (BMSPs) represent the advanced 
phase of stellar evolution in close, interacting binaries. Their 
observed orbital and stellar properties are fossil records of 
their evolutionary history. Thus one can use binary pul- 
sar systems as key probes of stellar astrophysics. It is 
well established that the neutron star component in bi- 
nary millisecond pulsar systems forms first, descending 
from the initially more massive of the two binary zero- 



age main sequence (ZAMS) stars. The neutron star is sub- 
sequently spun up to a high spin frequency via accre- 
tion of mas s and angular momentum once the secondary 
star evolves (lAlpar et al.[[ l982: Radhakr ishnan fc Srinivasanl 



1982; iBhattacharva fc van den Heuvellll991f ). In this recy 
cling phase the system is ob s ervable as a l ow-mass X-ray 
binar y (e.g. lHavakawal 1 19851 : iNagasd Il989l : iBildsten et"atl 
Il997h and t owards the end of this phas e as a n X-ray millisec- 
ond pulsar ijWi inands & van der Klislll998l : lArchibald et al.l 
l2009f: iPatruno fc Wattsll2012l h Although this standard for- 



© 0000 RAS 



2 Tauris, hanger & Kramer 



mation scenario is commonly accepted many aspects of the 
mass-transfer proces s and the accretion physics are still not 
understood in detail (|Lewin fc van der Klisl2006T ). Examples 
of such ambiguities include the accretion disk structure, the 
disk-magnetosphere transition zone, the accretion efficiency, 
the decay of the surface B-field of the neutron star and the 
outcome of common envelope evolution. 

The current literature on the recycling process of pul- 
sars is based on a somewhat simplified treatment of both the 
accretion process and the so-called equilibrium spin theory, 
as well as the application of a pure vacuum magnetic dipole 
model for estimating the surface B-field strength of a radio 
pulsar. These simplifications become a problem when trying 
to probe the formation and the evolution of observed recy- 
cled radio pulsars located near the classical spi n-up line for 
Eddington accretion in the PP-diagram (e.g. iFreire et al.l 
l201ll ). In this paper we discuss the concept and the location 
of the spin-up line and investigate to which extent it depends 
on the assumed details of the accretion disk/magnetosphere 
interactions and the magnetic inclination angle of the pulsar. 
Furthermore, we include the pla sma filled magnetosphere 
description by ISpitkovskvl J2006) for determining the sur- 
face B-field strengths of pulsars in application to spin-up 
theory and compare to the case of using the standard vac- 
uum dipole model. 

A key quantity to understand the formation of any given 
recycled radio pulsar is its spin period which is strongly re- 
lated to the amount of mass accreted. The amount of ac- 
cumulated mass is again determined by the mass-transfer 
timescale of the progenitor X-ray binary (and partly by the 
mode of the mass transfer and mass loss from the system), 
and hence it depends crucially on the nature of the donor 
star. At present the scenario is qualitatively well understood 
where low-mass X-ray binaries (LMXBs) and intermediate- 
mass X-ray binaries (IMXBs) in general lead to the for- 
mation of pulsars with a helium white dwarf (He WD) 
or a carbon-ox ygen white dwarf (CO WD) companion, re- 
spectively (e .i 



tigated in a future paper. For the evolution of ultra- 
compact X-ray binaries (UCXBs) leading to AXMSPs, and 
later to the formation of tight-orbit recycled radio BM- 
SP s with ultra-light ( < 0.08 Mr) companions, we refe r 
to iPodsiadlowski et all (l2002h: iDelove fc Bildstenl (l2003h : 



Tauris et al.| [2000; Podsiadlowsk i et al ]|2002l ; 
Tauris! 1201 ll ~ However, here we aim to quan- 



;ively (e .g. 
lDelovdl200S| ;rf 
tify this picture in more detail by re-analysing the spin-up 
process. As an example, we investigate in this paper the pos- 
sibilities of spinning up a neutron star from Case BB RLO 
leading to a mildly recycled pulsar with a CO WD compan- 
ion. We also discuss the increasing number of systems which 
apparently do not fit the characteristics of these two main 
populations of binary pulsars. 

We follow the pulsar spin evolution during the final 
phase of the mass transfer by including calculations of the 
torque acting on an accret ing pulsar du ring the Roche-lobe 
decoupling phase (RLDP; [Taurisll2012l ') for both an LMXB 
and an IMXB and compare the results. To complete the full 
history of pulsar spin evolution we subsequently consider the 
true ages of recycled radio pulsars (which are intimately re- 
lated to their spin evolution) by calculating isochrones and 
discussing the distribution of recycled pulsars in the PP- 
diagram. 

Although accreting X-ray millisecond pulsars (AXM- 
SPs) are believed to be progenitors of BMSPs, all of 
the 14 observed AXMSPs have orbital periods less than 
one day whereas fully recycled radio BMSPs are observed 
with orbital p eriods all th e way up to a few hundred 
days (see also iTaurij l2012t ). This puzzle will be inves- 



Ivan der Sluvs et al.l (|2005h ; Ivan Haaften et al.l (|2012h . In 
this paper we focus on the formation and evolution of BM- 
SPs with regular He WDs and CO/ON eMg WDs. 

The discovery of PSR J1614-2230 (|Hessels et al.ll2005l ; 
iDemorest et alj I2010T ) plays an important role for under- 
standing BMSP formation. This pulsar system is interest- 
ing since it has a unique combination of a short pulsar 
spin period of 3.2 ms and a massive CO WD companion. 
A rapidly spinning pulsar is usually associated with a long 
phase of mass-transfer evolution in an LMXB, whereas a 
CO WD companion in a relatively close orbit (8.7 days) is 
evidently the outcome of an IMXB evolution. A possibly 
solution to this paradox is that PSR J1614— 2230 evolved 
from Case A RLO of an IMXB which results in both a rel- 
atively long-lived X-ray phase (> 10 7 yr), needed to spin up 
the pulsar effectively, and leaving behind a CO WD. In 
this case PSR J1614-2230 is the fir st BMSP known t o 
have evolved via this path. Indeed, in iTauris et ail <|201ll ). 
hereafter Paper I, we investigated the progenitor evolu- 
tion of PSR J1614— 2230 with emphasis on the X-ray phase 
where the binary contained a neutron star and a donor star. 
We found two viable possibilities for the formation of the 
PSR J1614-2230 system: either it contained a 2.2 - 2.6 M© 
asymptotic giant branch donor star and evolved through a 
common envelope and spiral-in phase initiated by Case C 
RLO, or it descended from a close binary system with a 
4.0 — 5.0 Mq main sequence donor star via Case A RLO 
as hinted abov e. The latter scenario was also found by 
iLin et all (120111 ) . The fact that PSR J1614-2230 was spun- 
up to (less than) 3.2 ms could indeed hint which one of the 
two formation scenarios is most likely. In order to test this 
idea and to further distinguish between Case A and Case C 
we turn our attention, here in Paper II, to the spin dynamics 
of this pulsar in the recycling process. 

As discussed in Paper I, the distribution of neutron star 
birth masses is an important probe of both stellar evolu- 
tion, binary interactions and explosion physics. For a num- 
ber of BMSPs we are now able to put constraints on the 
birth mass of the pulsar given the derived amount of mass 
needed to spin up the observed recycled pulsar. In partic- 
ular, it is of interest to see if we can identify further pul- 
sars showing evidence of being born massive (~ 1.7 Mq) like 
PSR J1614-2230. 

In order to understand the many different observa- 
tional properties of BMSPs we have combined a variety of 
subtopics here with the aim of presenting a full picture of the 
subject. Given the many facets included, our new findings 
and the resulting length of the manuscript, we have chosen 
throughout this paper to finalize individual subtopics with 
a direct comparison to observational data followed by dis- 
cussions in view of our theoretical modelling. 

Our paper is structured in the follow way: We begin 
with a brief, updated review of the formation channels of 
BMSPs with CO WDs (Section [2} and present a summary 
of the latest observational data in Section In this section 
we also demonstrate an emerging unified picture of pulsar 
formation history which, however, is challenged by a num- 
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ber of interesting systems which share expected properties 
of both an LMXB and an IMXB origin. In Section 0] we 
investigate the recycling process in general with a focus on 
the location of the spin-up line in the PP-diagram and also 
relate the initial spin of a rejuvenated pulsar to the amount 
of mass accreted. The theoretical modelling is continued in 
Section [5] where we highlight the effects of the Roche-lobe 
decoupling phase on the spin evolution of recycled pulsars. 
In Section [B] our results are compared to the observed spin 
period distributions and in Section [7] we investigate if BM- 
SPs with CO WD companions obtained their fast spin pe- 
riods after a common envelope evolution. In Section [5] we 
discuss our results in a broader context in relation to the 
spin evolution and the true ages of millisecond radio pul- 
sars. In Section [5] we continue our discussion from Paper I 
on the formation and evolution of PSR J1614— 2230 and in 
Section [10] we return to the question of neutron star birth 
masses. Our conclusions are summarized in Section 1111 Fi- 
nally, in the Appendix we present a new tool to identify the 
most likely nature of the companion star to any observed 
binary pulsar. 



2 FORMATION OF BMSPS WITH CO WD 
COMPANIONS 

According to stellar evolution theory one should ex- 
pect radio pulsars to exist in binaries with a vari- 
ety of different companions: white dwarfs (He WDs, 
CO WDs and ONeMg WDs), neutron stars, black 
holes, ultra light semi-degenerate dwarfs (i.e. substel- 
lar companions or even planets), helium stars, main 
sequence stars an d, for very wide systems, su b -giant 
and giant stars dBhattacharva fc van den Heuvell 1 199 it 



iTauris fe van den Heuvell l200d . and references therein). Of 
these possibilities, pulsars orbiting black holes, helium stars 
and giant stars still remain to be detected. 
The majority of radio BMSPs have He WD companions. 
The formation of these systems is mainly channeled 
through LMXBs and have been well i n vestigated in previ- 
ous s t udies ( e.g. Webbink et "all 19831 ; Pvlvser fc Sayonfk] 



19881. 1 19891; iRappaport et al.l Il995l; lErgma et alj 



Tauris & Savonij el 1 19991: "IPodsiadlowski et al 



1998 



2002 



20041 ; Ivan der Sluvs et all 120051 : iDelovel 



Nelso n et al.l 

20081 ). These systems have orbital periods between less 
than 0.2 days and up to several hundred days. One of 
the most striking features of these systems is the relation 
between white dwarf masses and orbital periods (e.g. 
ITauris fc Savoniiel [19991 ). However, BMSPs with He WD 
companions may also form in IMXBs if Case A RLO is not 
initiated too late during the m a in sequence evolution of th e 
donor star (|Tauris et al.1 12000| ; IPodsiadlowski et all 120021 ). 
In Section [3] we identify 6 BMSPs which may have formed 
via this channel. 

BMSP systems with relatively heavy WDs (CO or 
ONeMg WDs) all have observed orbital periods ^ 40 days. 
For these systems there are a number of suggested formation 
channels which are briefly summarized below - see Paper I, 
and references therein, for a more rigorous discussion. In or- 
der to leave behind a CO WD, its progenitor must be more 
massive than 3.0 Mq if Roche-lobe overflow (RLO) is initi- 
ated while the donor star is still on the main sequence. If the 



donor star is an asymptotic giant branch (AGB) star when 
initiating mass transfer it can, for example, have a mass as 
low as 2.2 M e and still leave behind a 0.5 Mq CO WD rem- 
nant - see fig. 1 in Paper I. CO WDs in close binaries have 
masses between 0.33 — 1.0 Mq, whereas the ONeMg WDs are 
slightly heavier with masses of 1.1 — 1.3 Mq. The upper limit 
for the initial mass of donor stars in X-ray binaries leaving 
an ONeMg WD is often a ssumed to be about 7 — 8 Mq 
l|Podsiadlowski et al.l I20041 ). depending on the orbital pe- 
riod, metallicity, treatment of convecti on and the amount 
of con vective overshooting. However, IWellstein fc Langerl 
l|l999l) demonstrated that even stars with an initial mass in 
excess of 13 Mq may form an ONeMg WD if these progen- 
itor stars are in a tight binary and thus lose their envelope 
at an early stage via Case A RLO. Stars exceeding the up- 
per threshold mass for producing an ONeMg WD will leave 
behind a neutron star r emnant via an electron ca pture su- 
pernova (|Nomotol 1 19841 ; IPodsiadlowski et al.ll2004l ) or via a 
core-collapse supernova of type Ib/c for slightly more mas- 
sive stars. For single stars, or very wide binaries, the critical 
ZAMS mass for neutr on star production via type II super- 
novae is about 10 Mq jZhang et al.ll2008l ). 
In this paper we consider all systems with either a CO or an 
ONeMg WD as one population, denoted by CO WDs un- 
less specified otherwise, and all scenarios listed below apply 
to BMSPs with both types of massive white dwarf compan- 
ions. However, the required progenitor star masses are, in 
general, somewhat larger for the ONeMg WDs compared to 
the CO WDs. Despite these minor differences, the bottom 
line is that in all cases IMXBs are the progenitor systems of 
BMSPs with massive white dwarf companions. 



2.1 IMXBs with Case A RLO 

BMSPs with CO WDs, which formed via Case A RLO 
in IMXBs, evolve from donor stars w ith masses 3 — 5 Mq 
and orbital periods of a few days l|Tauris et al.l l|2000l) ; 
Podsiadlo wski et al.l l|2002l ); Paper I). The final orbital peri- 
ods of the BMSPs are 5 — 20 days. The mass-transfer phase 
consists of three parts: Al, A2 and AB. Phase Al is thermal 
time-scale mass transfer and lasts for about 1 Myr. The ma- 
jority of the donor envelope is transfered at a high rate ex- 
ceeding 10 -5 Mq yr _1 . However, the accretion onto the neu- 
tron star is limited by the Eddington luminosity, correspond- 
ing to an accretion rate of a few 10 -8 -Mq yr -1 , and thus 
the far majority of the transfered material (up to 99.9%) is 
ejected from the system. Phase A2 is driven by the nuclear 
burning of the residual hydrogen in the core. The result- 
ing mass-transfer rate is very low (< 10~ 9 M© yr -1 ) and 
only a few 10 -2 Mq of material is transfered towards the 
accreting neutron star in a time-interval of 20 — 50 Myr. 
Phase AB is caused by expansion of the donor star while un- 
dergoing hydrogen shell burning. The mass-transfer rate is 
10" 8 - 10" 7 Mq yr" 1 and lasts for ~ 10 Myr. A few 0.1 Mq is 
transfered towards the accreting neutron star - significantly 
more than in the previous two phases, Al and A2. As we 
shall see later on, this phase is responsible for spinning up 
the neutron star to become a millisecond pulsar. 
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2.2 IMXBs with early Case B RLO 

Wider orbit IMXB systems evolve from donor stars which 
are in the Hertzsprung gap at the onset of the RLO. The 
orbital periods at the onset of the Case B RLO are about 3 — 
10 days and donor masses of 2.5 — 5.0 Mq are needed to yield 
a CO WD remnant in a B MSP system (|Tauris et al.ll200d ; 
Podsiadlo wski et al.ll2002h . The final orbital periods of the 
BMSPs are about 3 — 50 days. The mass-transfer rate in the 
IMXB phase is highly super-Eddington (10 -5 Mq yr _1 ) and 
thus the far majority of the transfered material is ejected. 
However, a few 10 -2 Mq are transfered to the neutron star 
- enough to create a mildly (> 10 ms) spun- up millisecond 
pulsar, see Section!?] 

Binaries with either more massive donor stars (> 5 Mq) or 
initial periods in excess of 10 days will undergo dynamically 
unstable RLO (see Paper I, and references therein) leading 
to a common envelope and spiral-in. The outcome of the in- 
spiral is probably a merger for Case B RLO, since for these 
stars the (absolute) binding energy of the envelope is too 
large to allow for its ejection - except if accretion luminosity 
of the neutron star can be efficiently converted into kinetic 
energy of the envelope. 



2.3 IMXBs with Case C RLO and a common 
envelope 

Donor stars in systems with very wide orbits (Porb — 10 2 — 
10 3 days) prior to the mass-transfer phase develop a deep 
convective envelope as they become giant stars before filling 
their Roche-lobe. The response to mass loss for these stars 
is therefore expansion which causes the stars to overfill their 
Roche-lobe even more. To exacerbate this problem, bina- 
ries also shrink in size if mass transfer occurs from a donor 
star somewhat more massive than the accreting neutron 
star. This causes further overfilling of the donor star Roche- 
lobe resulting in enhanced mass loss etc. This situation is 
a vicious circle that leads to a dynamically unstable, run- 
away mass transfer and the formation of a common envelop e 
(CE) followed by a spiral-in phase , e.g. Paczvnskil (|l97rt ). 
Ilben fc Liviol l|l993h . Ilvanova et all (|2012T i. Whereas donor 
stars which initiate RLO during hydrogen shell burning still 
have relatively tightly bound envelopes, these stars here 
which in itiate RLO on the AGB have only wea kly bound en- 
velopes (|Han et al.lll994l ; lDewi fe Taurisil2000t ). Hence, these 
donor stars can survive the spiral-in phase of the neutron 
star without merging. To leave a CO WD remnant from 
such wide-orbit systems the ZAMS masses of the progenitor 
stars (the donor stars of the IMXBs) must be 2.2 — 6 Mq 
(Paper I), depending on the assumed amount of convective 
core-overshooting in the stellar models. 



2.4 BMSPs with a CO WD from LMXBs 

In addition to the formation channels from IMXBs, pulsars 
with a CO WD companion can, in rare cases, be produced 
from late Case B RLO in an LMX B system - e.g. see Ta- 
ble Al in lTauris fc Savoniid (1 19991 ). The outcome is an ex- 
tremely wide-orbit pulsar system (P or b > 1000 days) with a 
0.47 - 0.67 Mq CO WD and a slowly spinning pulsar. The 
pulsar B0820+02 (P = 0.86 sec, P orb = 1232 days) is an 



example of a system which followed this formation chan- 
nel (cf. Section f3.2|) . It is an interesting fact that CO WDs 
produced from donor stars in IMXBs can be less massive 
than CO WDs produced in such LMXBs. In low-mass stars 
(^ 2.2 Mq) helium is ignited in a flash under degener- 
ate conditions and CO WDs produced from these stars in 
LMXBs have minimum masses of 0.47 Mq, whereas in cores 
of intermediate-mass stars helium ignites non-degenerately 
and CO WDs from IMXBs can be made with masses down to 
0.33 Mm llKippenhahn fc Weigerd Il99d ; iTauris et all |2000| ; 
IPodsiadlowski et al.ll2002h . 

2.5 BMSPs with a CO WD from AIC? 

Under certain circumstances is may be possible for an accret- 
ing O-Ne-Mg WD to reach the Chandrasekhar mass limit 
and thereby implode to form a neutron star - the so-called 
accre t ion-induced collap se, AIC (|Mivaii et al.lll980l ; lNomotol 
1 19841 ; ICanal et al"1 [l990), A neutron star formed this way 
may leave behind a radio pulsar. The necessary conditions 
for a pulsar formed this way to further accrete from the 
companion star and leaving behind, for example a BMSP 
with a CO WD companion, is investigated in a separate pa- 
per (Tauris et al., in preparation). However, the recycling 
process, which is the main topic in this paper, remains the 
same. 



3 OBSERVATIONAL CHARACTERISTICS OF 
BMSPS WITH CO WD COMPANIONS 

Before looking into the details of recycling a pulsar in an 
X-ray binary, we first summarize the observational proper- 
ties of the end products, the BMSPs. This is important for 
our subsequent discussions and understanding of the the- 
oretical modelling which constitutes the core of our work. 
Furthermore, we shall investigate if and how the BMSPs 
with CO WD companions can be distinguished from the 
large population of BMSPs with He WD companions. For 
this purpose we make use of the PP-diagram, the Corbet di- 
agram and information of orbital eccentricities. 

3.1 The PP-diagram 

In order to understand the formation and evolution of recy- 
cled pulsars we have plotted in the PP-diagram in Fig. [T] 
the distribution of 103 binary radio pulsars located in the 
Galactic disk. The various companion types are marked 
with different symbols (see also Table [1} ■ The majority of 
the companion stars are helium white dwarfs (He WDs). 
However, there is also a significant population of the more 
massive CO and ONeMg WDs, as well as pulsars with 
another neutron star as companion. In some systems the 
companion is still a main sequence star. In other (very 
tight) systems the companion has evolved into a semi- 
degenerate substellar remnant (< 0.08 Mq) - possibly aided 
by evaporation caused by the illumination from the pul- 
sar fo l lowing an ultra-compact LMXB phase (IPod siadlowski 
199ll; iDeloye fc Bildstenl l2003l; Ivan der Sluvs et all 120051 ; 



Bailes et alj|201ll ; Ivan Haaften et al.ll2012r i. 
Pulsar systems found in globular clusters are not useful as 
probes of stellar evolution and binary interactions since their 
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Spin period, P (sec) 

Figure 1. Distribution in the PP-diagram of 103 binary radio pulsars in the Galactic disk. The five solid lines indicate theoretical spin- 
up lines using equation dlOI I for a 1.4 Mq pulsar accreting with a rate of M = M^ddi 0.1 -^Eddi 0-01 AfEddi 0.001 M and M = 10 -4 MBddj 
respectively (top to bottom). Also shown as dotted lines are characteristic ages of 100 Myr, 1 Gyr and 10 Gyr, respectively. The location 
of PSR J1614— 2230 is peculiar in the sense that it has a massive CO WD companion (blue diamond) in a region where all other pulsars 
have a low-mass He WD companion (red filled circles). The six pulsars marked with circles may have evolved from IMXBs and some 
of these are potential candidates for having a CO WD companion — see Section 13.21 Data taken from the ATNF Pulsar Catalogue, 
|http: / /www.atnf.csiro.au/research/pulsar/psrcat/l ^Manchester et al,ll2005T l. in May 2012. 



location in a dense environment causes frequent exchange 
collisions whereby the orbital parameters are perturbed or 
the companion star is replaced with another star and thus 
information is lost regarding the mass transfer and spin-up 
of the pulsar. These pulsars could even have undergone ac- 
cretion events from two or more different companions. 

The location of PSR J1614-2230 in the PP-diagram is 
quite unusual for a BMSP with a CO WD companion. It is 
located in an area which is otherwise only populated by BM- 
SPs with He WD companions. This is important as it proves 
for the first time that efficient, full recycling (i.e. P ^ 8 ms) 
of a pulsar is also possible in X-ray systems leading to BM- 
SPs with CO WD companions. In general, the BMSPs with 
CO WD companions are seen to have large values of both P 
and P compared to BMSPs with He WD companions. This 
trend can be understood since the former systems evolve 
from IMXBs which often have short lasting RLO and thus 
inefficient spin-up - PSR J1614— 2230 being the only known 
exception (this system evolved from an IMXB via Case A 
RLO, see Paper I). 

In Table [2] we list all binary radio pulsars with CO WD 
(or ONeMg WD) companions - see also lHobbs et al l l|2004h . 
At the bottom of the Table we list an additional number of 
sources which we identify as CO WD candidate systems (al- 
ternatively, these candidates may be He WDs evolving from 



Table 1. Binary radio pulsars and their companions. Data taken 
from the ATNF Pulsar Catalogue in May 2012. 



Population Number 



All known binary pulsars: 186 

Galactic disk 110 

Extra galactic (SMC) 1 

Globular clusters 75 



Galactic disk* binary pulsars with measured P: 103 

Companion star: 

He white dwarf 56 

CO/ONeMg white dwarf 19 

Neutron star 10 

Main sequence star* 6 

Ultra light (semi)deg. dwarf, or planet(s) 12 



* Including PSR J0045-7319 which is locate d in the Sma ll 
Magellanic Cloud and has a B1V companion jKaspi et alJll994h . 



IMXBs). The three slowest rotating of the 20 pulsars with a 
CO WD companion are not recycled - either becau se they 
were formed after the WD |Tauris fc Sennelsl :2000). or be- 
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Table 2. Binary radio pulsars which are likely to have evolved 
from IMXB systems. Data taken from the ATNF Pulsar Cata- 
logue in May 2012. 



PSR 




r orb 


p 


p 




JW WD 






days 




1Q-19 


10 -4 


Mr-, 
1V1Q 


J1952- 


-2630 


0.392 


20.7 


n/a 


n/a 


1.13 


J1757- 


-5322 


0.453 


8.87 


0.263 


0.0402 


0.67 


J1802- 


-2124 


0.699 


12.6 


0.726 


0.0247 


0.78* 


B0655- 


h64 


1.03 


196 


6.85 


0.0750 


0.80 


J1435- 


-6100 


1.35 


9.35 


0.245 


0.105 


1.08 


J1439- 


-5501 


2.12 


28.6 


1.42 


0.499 


1.30* 


J1528- 


-3146 


3.18 


60.8 


2.49 


2.13 


1.15 


J1157- 


-5112 


3.51 


43.6 


1.43 


4.02 


1.30* 


J1337- 


-6423 


4.79 


9.42 


1.95 


0.197 


0.95 


J1603- 


-7202 


6.31 


14.8 


0.156 


0.0928 


0.34 


J2145- 


-0750 


6.84 


16.1 


0.298 


0.193 


0.50 


J1022-f 


-1001 


7.81 


16.5 


0.433 


0.970 


0.85 


J062H 


-1002 


8.32 


28.9 


0.473 


24.6 


0.67* 


J1614- 


-2230 


8.69 


3.15 


0.0962 


0.0130 


0.50* 


J1454- 


-5846 


12.4 


45.2 


8.17 


19.0 


1.05 


J0900- 


-3144 


18.7 


11.1 


0.491 


0.0103 


0.42 


J1420- 


-5625 


40.3 


34.1 


0.675 


35.0 


0.44 


J1141- 


-6545° 


0.198 


394 


4307 


1719 


1.02* 


B2303- 


K6 a 


12.3 


1066 


569 


6584 


1.30* 


B0820- 


402 b 


1232 


865 


105 


118.7 


> 0.52 


J1622- 


-6617 c 


1.64 


23.6 


0.636 


n/a 


0.11 


J1232- 


-6501 c 


1.86 


88.3 


8.11 


1.09 


0.17 


J1745- 


-0952 c 


4.94 


19.4 


0.925 


0.0985 


0.13 


J18414 


-0130 c 


10.5 


29.8 


81.7 


0.819 


0.11 


J19044 


-0412 c 


14.9 


71.1 


1.10 


2.20 


0.26 


J1810- 


-2005 c 


15.0 


32.8 


1.47 


0.192 


0.33 



* Mass obtained from timing. Other masses are median masses 
(i.e. calculated for i = 60° and an assumed A/ns = 1.35 Mq). 
a A non-recycled pulsar formed after its WD companion 
iTauris fc SenneIsj|20o"oT) . 

b This system descends from a very wide-orbit LMXB (see text). 
c CO WD candidate or He WD evolving from an IMXB. 

cause they formed in a very wide orbit (cf. Section f2 . 4 [) . Of 
the remaining 17 recycled pulsars with a CO WD compan- 
ion, 16 have a measured value of P. These are plotted in 
Fig. [2] together with those other systems which may also 
have evolved from IMXB systems and thus possibly host a 
CO WD companion too (cf. Section [3.2[) . With the excep- 
tion of PSR J1841+0130 (which is young, see Section lSXI)) . 
these pulsars are distributed in a somewhat more linear man- 
ner in the PP-diagram compared to BMSPs with a He WD 
companion. 



3.2 The Corbet diagram - six unusual pulsars 

Plotting BMSPs in the Corbet diagram (|CorbetHl984T ) could 
lead to potential clues about their formation history. In the 
top panel of Fig. we have shown the distribution of all 
binary pulsars. Little can be learnt from this plot when con- 
sidering the entire ensemble of binary pulsars as one pop- 
ulation. In the central plot, however, we display only those 
binary pulsars which have a He WD companion. The dis- 
tribution of pulsars in this diagram shows some interesting 
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Figure 2. The observed recycled pulsars with CO WD compan- 
ions (diamonds) seem to cluster somewhat along a straight line 
in the PP— diagram. The full line is a linear regression fit to all 
these pulsars and the dotted lines represent constant character- 
istic ages. Other pulsars which are also candidates for having 
an IMXB origin are shown as dots in a circle (cf . Section 13-21 ) . 
PSR J1841+0130 is a young example of such a candidate — see 
text. 



features. The six pulsars marked with a circle are noticeable 
for having unusual slow spin periods compared to other pul- 
sars with similar orbital periods, see Table [5] for further de- 
tails. This fact hints an inefficient and sh orter sp i n-up phase 
compared to the BMSPS with fast spin. iTaurisI (l201ll ) sug- 
gested that these pulsars may perhaps have formed via the 
accretion induced collapse (AIC) of a WD. Here, we suggest 
instead that these pulsars (which all have P or b 1 day) 
may have formed in IMXBs. In that case the companions are 
likely to be CO WD s, although He WD companions may also 
form from IMXBs (|Tauris et al.ll2000| ; iPodsiadlowski et all 
2002) . Similar ideas of an IMX B orig in have been proposed 
bv lCamilo et all (|200ll ) andQ (|2002l ). The latter author ar- 
gued that accretion disk instabilities in IMXBs may explain 
the slow spin periods and the high P values of these pulsars. 
The pulsars with P or b > 200 days all have slow spin peri- 
ods - see Section T6, 31 for a discussion on these systems. The 
bottom panel shows the binary pulsars with CO WD com- 
panions. The main thing to notice is that these systems have 
Porb ^ 40 d and slower pulsar spin periods compared to BM- 
SPs w ith He WD companions, as first pointed out bv lCamilol 
(1996). The slower spin of pulsars in this population can be 
understood well which is discussed further in Section [6] A 
few pulsars marked with a circle or a triangle stand out 
from the rest of the population. These are PSR J1141— 6545 
and PSR B2303+46 whe re the neutron star form ed after the 
white dwarf companion (jTauris fc Senneldl200ch and there- 
fore these pulsars remain non-recycled with slow spin peri- 
ods, and PSR B0820+02 which formed from an extremely 
wide-orbit LMXB with short lasting and possibly inefficient 
RLO (see Sections [21] and [OJ. 
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Figure 3. All binary radio pulsars in the Galactic disk plotted 
in a Corbet diagram (top), see Fig. [JJ for an explanation of the 
symbols plotted. Due to their missing P measurements 8 of these 
pulsars were not shown in Fig. [T] The central and bottom panels 
show the distribution of pulsars with He WD (60 systems) and 
CO WD companions (20 systems), respectively. In the bottom 
panel the grey shaded region marks the range of spin periods 
obtained from one of our model calculations of Case BB RLO in 
a post-CE binary (Fig. [lOj . See text for further discussion. 
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Figure 4. Eccentricity as a function of orbital period for bi- 
nary radio pulsars in the Galactic disk. The various symbols are 
defined in Fig. [TJ We notice PSR J1822-0848 which has an ec- 
centricity two orders of magnitude larger than other wide-orbit 
binary pulsars with a He WD companion. 



3.2.1 PSR J1841+0130 

PSR J1841+0130 (discovered by lLorimer et alll2006] ) is one 
of the six pulsars which we note to have an offset location in 
the Corbet diagram when comparing to BMSPs with He WD 
companions (Fig. [3] central panel) . It has a slow spin period 
of 29 ms and an orbital period, P or b = 10. 5 d . The majority 
of other BMSPs with similar P or b near 10 days have much 
faster spin periods between 2—10 ms. The slow spin pe- 
riod of PSR J1841+0130 can be explained if it originated 
from an IMXB system rather than an LMXB system. The 
mass function of this pulsar, / = 4.22 x 10 -4 Mq suggests 
that its companion star has a mass of 0.11 Mq for an or- 
bital inclination angle of 60°. If the companion star is a 
CO WD formed in an IMXB it must have a mass ^ 0.33 Mq 
(cf. Section T2.4J1 . This would require a small inclination an- 
gle, i ^ 20° (depending on the neutron star mass). The 
probability for this to happen from a distribution of random 
orbital orientations is about 6 per cent. The situation is sim- 
ilar in PSR J1622— 6617 which is another candidate systems 
to have evolved from an IMXB. However, IMXBs can also 
leave behind He WDs which have somewhat smaller masses. 
Alternatively, if PSR J1841+0130 originated from an LMXB 
system then according to the (Mm), P rb)-relation, e.g. 
iTauris fc Savoniiel (|l999r i. we would expect i w 25° to obtain 
the predicted He WD mass of 0.26 ± 0.01 Mq. (The proba- 
bility for this is about 9 per cent.) Hence, even in this case 
the inclination angle would be rather low. As discussed in 
Sections 1431 and lOH PSR J1841+0130 is also interesting 
for its ability to constrain spin-up physics and for its young 
age. 



3.3 The orbital eccentricity 

Anothe r fossil record of binary evolution is the orb ital eccen- 
tricity (|PhinneJl99l iPhinnev fc Kulkarnilll994l ). In Fig.H 
we have plotted the eccentricity as a function of orbital pe- 
riod. The BMSPs with He WDs have in general somewhat 
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lower eccentricities than systems with CO WDs. The spread 
is large for both of the two populations which also over- 
lap. Among the systems with CO WD companions it is not 
surprising that PSR J1614— 2230 has the lowest eccentric- 
ity since we have demonstrated (Paper I) that this system 
formed via a relatively long phase of stable RLO. It is in- 
teresting to notice that the four BMSPs with CO WDs and 
Porb < 2 days all have quite small eccentricities ^ 10 even 
though they are believed to have formed via a CE. 



3.3.1 PSR J '1822- 0848 

PSR J1822-0848 (lLorimer et al.ll2006l ) has an orbital period 
of 287 days and an eccentricity, ecc = 0.059. From Fig. 2] 
it is seen that this eccentricity is two orders of magnitude 
larger than that of other pulsars in similar wide-orbit sys- 
tems. This fact, combined with its very slow spin period of 
835 ms (see the central panel in Fig. [3} and its high value 
of P = 1.35 x 10~ 16 , hints that this pulsar may not have 
a He WD companion star. It is possible it belongs to the 
same class as PSR B0820+02 which is in a wide orbit with 
a CO WD (see Section [2. 4|l . However, in that case it needs 
to be explained why some pulsars have He WD companions 
in even larger orbits compared to PSR J1822— 0848. Alter- 
natively, one could speculate that PSR J1822— 0848 experi- 
enced weak spiral-in from an almost unbound common en- 
velope of an AGB star. This could explain its non-recycled 
characteristics and why it has P or b < 1000 days. 
Although, it has a relatively large eccentricity of 0.059 this 
value is still far too small to suggest that PSR J1822— 0848 
was born after the W D, like in the case of PS R 1141—6545 
and PSR B2303+46 (|Tauris fc Sennelsl 120001 ). Even in the 
unrealistic case of a collapse of a naked core with no enve- 
lope, and without a kick, the present eccentricity of 0.059 
for such a wide orbit would require a maximum instant mass 
loss of only ~ 0.10 Mq - a value which is even less than the 
release of the gravitational binding energy during the core 
collapse. Thus this scenario is not possible and we therefore 
conclude that PSR J1822-0848 was formed before its WD 
companion, like 99 per cent of all detected binary pulsars 
with WD companions. 



4 RECYCLING THE PULSAR 

We now proceed with a general examination of the recycling 
process of pulsars. We begin by analysing the concept of 
spin-up lines in the PP-diagram and derive an analytic 
expression for the amount of accreted mass needed to 
spin up a pulsar to a given equilibrium spin period. We use 
standard theory for deriving these expressions, but include 
it here to present a coherent and detailed deriva tion needed 
for ou r purposes. Furthermore, we apply the S pitkovskvl 
(2006) solution to the pulsar spin-down torque which 
combines the effect of a plasma current in the magneto- 
sphere with the magnetic dipole model. The physics of 
the disk-magnetosphere interaction is still not known in 
detail. The interplay between the neutron star magnetic 
field and the conducting plasma in the accretion disk is a 
rather complex proces s. For further de t ails o f t he accretion 
iPringle fc Reesl d 19721); lLamb et al 



l|l979al lb h: IShapiro fc Teukolskvl dl983l); iGhosh fc LambI 
1 19921) nSpruit fc Taaml dl993h: ICampana et al J il 19981) ; 



Frank et all (I2002T): iRappaport et all (120041); iBozzo et al 



physic s 



we refer to ^_ 
iDavidson fc Ostrikerl 



( 1973f ); IGhosh fc L amb 



|2009i) ; TLVAngelo fc Spruitl l|2010l ); llkhsanov fc Beskrovnava 
(|2012i ) and references therein. 

4.1 The accretion disk 

The accreted gas from a binary companion possess large 
specific angular momentum. For this reason the flow of gas 
onto the neutron star is not spherical, but leads to the for- 
mation of an accretion disk where excess angular momen- 
tum is transport ed outwards by (turbulent-e n hanced) vis- 
cous stresses, e.g. IShapiro fc Teukolskvl (Il983l ); iFrank et al.1 
(2002). Depending on the mass-transfer rate, the opacity of 
the accreted material and the temperature of the disk, the 
geometric shape and flow of the material may take a va- 
riety of forms (thick disk, thin thick, slim disk, torus-like, 
ADA F). Popular models of the inner disk l|Ghosh fc LambI 
Il992l) include optically thin/thick disks which can be either 
gas (GPD) or radiation pressure dominated (RPD). The ex- 
act expression for the spin-up line in the PP-diagram also 
depends on the assumed model for the inner disk - mainly 
as a result of the magnetosphere boundary which depends 
on the characteristics of the inner disk. Close to the neu- 
tron star surface the magnetic field is strong enough that 
the magnetic stresses truncate the Keplerian disk, and the 
plasma is channeled along field lines to accrete on to the 
surface of the neutron star. At the inner edge of the disk 
the magnetic field interacts directly with the disk material 
over some finite region. The physics of this transition zone 
from Keplerian disk to magnetospheric flow is important 
and determines the angular momentum exchange from the 
differential rotation between the disk and the neutron star. 
Interestingly enough, it seems to be the case that the resul- 
tant accretion torque, acting on the neutron star, calculated 
using detailed models of the disk-magnetosphere interaction 
does not deviate much from simple expressions assuming 
idealized, spherical accretion and ne wtonian dynamics. For 
example, it has been pointed out bv lGhosh fc LambI (|l992T l 
as a fortuitous coincidence that the equilibrium spin period 
calculated under simple assumptions of spherical flow re- 
sembles the more detailed models of an optically thick, gas 
pressure dominated inner accretion disk. Hence, as a start- 
ing point this allows us to expand on standard prescriptions 
in the literature with the aims of: 1) performing a more 
careful analysis of the concept of a spin-up line, 2) deriv- 
ing an analytic expression for the mass needed to spin up 
a given observed millisecond pulsar, and 3) understanding 
the effects on the spinning neutron star when the donor star 
decouples from its Roche-lobe. As we shall discuss further 
in Section [5] the latter issue depends on the location of the 
magnetospheric boundary (or the inner edge of the accretion 
disk) relative to the corotation radius and the light-cylinder 
radius of the pulsar. All stellar parameters listed below will 
refer to the neutron star unless explicitly stated otherwise. 



4.2 The accretion torque — the basics 

The mass transfered from the donor star carries with it angu- 
lar momentum which eventually spins up the rotating neu- 
tron star once its surface magnetic flux density, B is low 
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enough to allow for efficient accretion, i.e. following initial 
phases where either the magneto-dipole radiation pressure 
dominates or propeller effects are at work. 
The accretion torque acting on the spinning neutron star 
has a contribution from both material stress (dominant 
term), magnetic stress and viscous stress, and is given by: 
N = J* = (d/dt)(ISl+), where J* is the pulsar spin an- 
gular momentum, fi* is its angular velocity and I ~ 1- 
2 x 10 45 g cm 2 is its moment of inertia. The exchange of an- 
gular momentum (J = fx p) at the magnetospheric bound- 
ary eventually leads to a gain of neutron star spin angular 
momentum which can approximately be expressed as: 



The result is: 



A J* = VGMr A AM £ 



(1) 



where £ ~ 1 is a numerical factor whic h depends on the flow 
pattern ijGhosh fc Lamblll979bl . 1 1992( 1. AM = M ■ At is the 
amount of mass accreted in a time interval At with average 
mass accretion rate M and 



''a 



/ B 2 R 6 \ 
\MV2Gm) 



2/7 



(2) 



22 km B 



4/7 



M 



0.1 M E 



-2/7 



M 



1.4 M© 



-5/7 



is the Alfven radius defined as the location where the mag- 
netic energy density will begin to control the flow of mat- 
ter (i.e. where the incoming material couples to the mag- 
netic field lines and co-rotate with the neutron star mag- 
netosphere). Here B is the surface magnetic flux density, 
R is the neutron star radius, M is the neutron star mass 
(see relation between R and M further below) and Bg is 
B in units of 10 s Gauss. A typical value for the Alfven ra- 
dius in accreting X-ray millisecond pulsars (AXMSPs), ob- 
tained from B ~ 10 s G and M ~ 0.01 M Ed d, is ~ 40km 
corresponding to about 3 R. The expression above is found 
by equating the magnetic energy density (B 2 /8iv) to the 
ram pressure of the incoming matter and using the conti- 
nuity equation (e.g. IPringle fc Reeslll972l ). Furthermore, it 
assumes a scaling with distance, r of the far-field strength 
of the dipole magnetic moment, \i as: B(r) oc fJ-/r 3 (i.e. dis- 
regarding poorly known effects such as magnetic screening 
l|Vasvliunaslll979D '). A more detailed estimation of the loca- 
tion of the inner edge of the disk, i.e. the coupling radius 
or magnetospheric bou ndary, is given by: r maK = if> ■ ta , 
where <t> is 0.5 - 1.4 (iGhosh fc Lambl 1 19921 : IWand 1 19971 : 



iD'Angelo fc Spruitll201ol ) 



4-2.1 The surface B-field strength of recycled radio pulsars 

Before we proceed we need an expression for B. One can 
estimate the B-field of recycled MSPs based on their ob- 
served spin period, P and its time derivative P. The usual 
assumption is to apply the vacuum magnetic dipole model 
in which the rate of rotational energy loss (E Iot = IQCl) is 
equated to the energy-loss rate caused by emission of dipole 
waves (with a frequency equal to the spin frequency of the 
pulsar) due to an inclined axis of the magnetic dipole field 
with respect to the rotation axis of the pulsar: 



dipole = (-2/3c 3 )|/i| 2 



B di 



pole 



3c 3 IPP 1 



8Tf 2 i? 6 sin a 



1.6 x 10 19 G 



(4) 



M 



1.4 M, 



3/2 



10 45 gem 2 , both of 



where the magnetic inclination angle is < a ^ 90° . This is 
the standard equation for evaluating the B-field of a radio 
pulsar. Our numerical scaling constant differs by a factor of 

3.2 x 10 1 



a few from the conventional one: B 
which assumes R — 10 km and I 
which we believe are slightly underestimated values (our as- 
sumed relations between 7, M and R are discussed in Sec- 
tion ^. 2. 2"j) . Note also that some descriptions in the literature 
apply the polar B-field strength (B p = 2 B) rather than the 
equatorial B-field strength used here. 

It is important to realize that the above expression does 
not include the rotational energy loss obtained when con- 
sidering the spin-down torque caused by the j x B force 
exerted by the plasma current in the ma gnetosphere (e.g. 
iGoldreich fc Julian 1969; Spitkovskv 2008). For this reason 
the vacuum magnetic dipole model does not predict any 
spin-down torque for an aligned rotator (a = 0°), which 
is not correct. The incompleteness of the vacuum magnetic 
dipole model was in part icular evident after the discovery 
of intermittent pulsars bv lKramer et al.1 (|2006l ) and demon- 
strated the need for including the plasma term in the spin- 
down torque. A combined model was derived bv lSpitkovskvl 
(2006) and applying his relation between B and a we can 
rewrite the above expression slightly: 



B 



/ c 3/pp 


1 


4ty 2 R 6 


1 + sin 2 a 



(5) 



M 



3/2 



1.4 M J V l + sin 2 Q 
This new expression leads to smaller estimated values of B 
by a factor of at least y/3, or more precisely y§(2~+ cot 2 a), 
compared to the vacuum dipole model. As we shall shortly 
demonstrate, this difference in dependence on a is quite 
important for the location of the spin-up line in the PP— 
diagram. 

4-2.2 The Eddington accretion limit 
The Eddington accretion limit is given by: 



Mi 



(7T 



^R^e (6) 



3.0 x 10" 8 M yr _1 JJi 



1.3 

1 + X 



where c is the speed of light in vacuum, 



is the pro- 
ton mass, <tt is the Thomson scattering cross section, /i c = 
2/(1 + X) is the mean molecular weight per electron which 
depends on the hydrogen mass fraction, X of the accreted 
material, and R\z is the neutron star radius in units of 
13 km. The expression is found by equating the outward ra- 
diation pressure to the gravitational force per unit area act- 
ing on the nucleons of the accreted plasma. In general, lumi- 
(3) nosity is generated from both nuclear burning at the neutron 
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star surface as well as from the release of gravitational bind- 
ing energy of the accreted material, i.e. L — (e nuc + e acc ) M. 
However, for accreting neutron stars e nuc <g e aC c and thus we 
have neglected the contribution from nuclear processing. To 
estimate the neutron star radius we used a mass-radius ex- 
ponent following a simple non-relativistic degenerate Fermi- 
gas polytrope (R oc Af~ 1//3 ) with a scaling factor such that: 
R = 15 (M /M a )~ 1/3 km, ca l ibrate d from PSR J1614-2230, 
cf. fig. 3 in lDemorest et al.l l|201Ch . In our code the value of 
MEdd is calculated according to the chemical composition of 
the transfered material at any given time. Note, the value of 
M^dd is only a rough measure since the derivation assumes 
spherical symmetry, steady state accretion, Thompson scat- 
tering opacity and Newtonian gravity. 



4.3 The spin-up line 

The observed spin evolution of accreting neutron stars of- 
ten shows rather sto chastic variations on a short timescale 
jBildsten et al.ll 19971 ). The reason for the involved dramatic 
torque reversals is not well known - see hypotheses listed 
at the beginning of Section [5] However, the long-term spin 
rate will eventually tend towards the equilibrium spin period, 
P Gq - meaning that the pulsar spins at the same rate as the 
particles forced to corotate with the B-field at the magne- 
tospheric boundary. The location of the associated so-called 
spin-up line for the rejuvenated pulsar in the PP-diagram 
can be found by considering the equilibrium configuration 
when the angular velocity of the neutron star is equal to the 
Keplerian angular velocity of matter at the magnetospheric 
boundary where the accreted matter enters the magneto- 
sphere, i.e. £1+ — Q ocl = WcSiK^mag) or: 



p 

i cq — 



2tt 



GM Lj c 



(7) 



1.40 ms Bt 



M 



-3/7 



AI 



1.4M P 



-5/7 



.O.lMEdd 

where 0.25 < uj c ^ 1 is the so-called critical fastness param- 
eter which is a measure of when the accretion torque van- 
ishes (depending on the dynamical import ance of the pul- 
sar spi n rate and the magnetic pitch angle, I Ghosh fc Lambl 
ll979bT ). One must bear in mind that factors which may 
differ from unity were omitted in the numerical expression 
above. In all numerical expressions in this paper we assumed 
sin a = cf> = £ = (jJc = 1- Actually, the dependence on the 
neutron star radius, R disappears in the full analytic for- 
mula obtained by inserting equations (J2]| and ([5]) into the 
top expression in equation l[7). and using r mag = ■ ta, 
which yields: 



R 18/7 



I (1 + sin a) 



-3/4 ,21/8 



-7/4 



(8) 
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cq \V2G 5 M 5 M 3 J 

Notice, in the above step we needed to link the B-fields of 
accreting neutron stars to the B-fields (expressed by P and 
P) estimated for observed recycled radio pulsars. This con- 
nection can be approximated in the following manner: If the 
radio pulsar after the recycling phase is "born" with a spin 
period, Po which is somewhat close to P cq (i.e. if the Roche- 
lobe decoupling phase (RLDP) did not significantly affect 
the spin period of the pulsar, cf. discussion in Section [5]), 
then we can estimate the location of its magnetosphere when 



the source was an AXMSP just prior to the accretion turn- 
off during the RLDP. (Further details of our assumptions 
of the B-fields of accreting neutron stars are given in Sec- 
tion 14.4.11 and in Section [8] we discuss the subsequent spin 
evolution of recycled radio pulsars towards larger periods, 
P > Po.) 

It is often useful to express the time derivative of the 
spin period as a function the equilibrium spin period, for 
example for the purpose of drawing the spin-up line in the 
PP-diagram: 



2 l/6 G 5/3 jVf M 5/3p4/3 

P = ~^JW^ J ( X + sm Q ) 



-7/2 7/3 



(9) 



Given that A^Edd is a function of the neutron star radius and 
using the relation between M and R stated in Section [47272] 
we need a relation between the moment of inertia and the 
mass of the neutron star. Accord ing to the equations-of-state 
studied by IWorlev et al.l (|2008h these quantities scale very 
close to linearly as: Z45 ~ M/Mq (see their fig. 4) where 
I45 is the moment of inertia in units of 10 45 g cm 2 . Towards 
the end of the mass-transfer phase the amount of hydrogen 
in the transfered matter is usually quite small (X ^ 0.20). 
The donor star left behind is basically a naked helium core 
(the proto WD). Hence, we can use equation ([6} and rewrite 
equation Q to estimate the location of the spin-up line for 
a recycled pulsar in the PP-diagram only as a function of 
its mass and its mass-accretion rate: 



P = 3.7 x 10" 



(M/M ) 



2/3 p 4/3 



M 

A/Ed, 



(10) 



assuming again sin a = cj> = ui c — 1, and where P ms is the 
equilibrium spin period in units of milliseconds. 

In the literature the spin-up line is almost always plot- 
ted without uncertainties. Furthermore, one should keep in 
mind the possible effects of the applied accretion disk model 
on the location of the spin-up line, cf. Section \4. II In Fig. [5] 
we have plotted equation <(9j for different values of a, <f> and 
lj c to illustrate the uncertainties in the applied accretion 
physics to locate the spin-up line. The upper boundary of 
each band (or "line") is calculated for a neutron star mass 
M = 2.0 Mq and a magnetic inclination angle, a = 90°. The 
lower boundary is calculated for M = 1.0 Mq and a = 0°. 
The green hatched band corresponds to <j> = 1 and uj c = 1, 
which we used in our calculations throughout this paper. 
The blue and red hatched bands are upper and lower limits 
set by reasonable choices of the two parameters (0, w c ). In 
all cases we assumed a fixed accretion rate of M = MEdd- 
The location of the spin-up line is simply shifted one order 
of magnitude in P down (up) for every order of magnitude 
M is decreased (increased). 

It is important to realize that there is no universal spin- 
up line in the PP-diagram. Only an upper limit. Any in- 
dividual pulsar has its own spin- up line/location which in 
particular depends on its unknown accretion history (M). 
Also notice that the dependence on the magnetic inclina- 
tion angle, a is much less pronounced when applying the 
Spitkovsky formalism for estimating the B-field of the MSP 
compared to applying the vacuum dipole model. The dif- 
ference in the location of spin-up lines using a = 90° and 
ct = 0° is only a factor of two in the Spitkovsky formalism. 
For a comparison, using the vacuum dipole model with a 
small magnetic inclination of a — 10° results in a spin-up 
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1 10 100 

Spin period, P (ms) 

Figure 5. The spin-up line is shown as three coloured bands 
depending on the parameters (<f>, u c ). In all three cases the 
spin-up line is calculated assuming accretion at the Eddington 
limit, M = MEdd and applying the Spitkovsky torque formal- 
ism. Thus these "lines" represent upper limits for the given set of 
parameters. The width of each line (band) results from using a 
spread in neutron star mass and magnetic inclination angle from 
(2.0 Mq, a = 90°) to (1.0 M©, a = 0°), upper and lower bound- 
ary, respectively. The dashed line within the upper blue band is 
calculated from (2.0 JW©, a = 0°). Hence, the band width above 
this line reflects the dependence on a, whereas the band width 
below this line shows the dependence on M. The two red dashed 
spin- up lines below the red band are calculated for a 1.4 M© neu- 
tron star using the vacuum magnetic dipole model for the radio 
pulsar torque and assuming cf> = 1.4, u) c = 0.25 and M = A^Edd 
for a = 90° (upper) and a = 10° (lower). The observed distribu- 
tion of binary and isolated radio pulsars in the Galactic disk (i.e. 
outside globular clusters) are plotted as filled and open circles, 
respectively - see Fig. 1 131 for further information. Also plotted is 
the pulsar J1823— 3021A which is located in the globular cluster 
NGC 6624. This pulsar and the three other pulsars marked by a 
red circle arc discussed in the text. 



line which is translated downwards by almost two orders 
of magnitude compared to its equivalent orthogonal rotator 
model, cf. the two red dashed lines in Fig. [5] 

If we assume that accretion onto the neutron star is 
indeed Eddington limited then the three bands in Fig. [S] 
represent upper limits for the spin-up line for the given sets 
of (</>, w c )- Thus we can in principle use this plot to constrain 
(4>, cj c ) and hence the physics of disk-magnetosphere interac- 
tions. The fully recycled pulsars B1937+21 and J0218+4232, 
and the mildly recycled pulsar PSR J1841+0130, are inter- 
esting since they are located somewhat in the vicinity of the 
green spin-up line. Any pulsar above the green line would 
imply that <f> < 1. The pulsar PSR J1823-3021A is close to 
this limit. Usually the derived value of P for globular clus- 
ter pulsars is influenced by the cluster potential. However, 
the P value for PSR J1823— 3021A w as recently constra ined 
from Fermi LAT 7-ray observations |Freire et al.ll201lh and 
for this reason we have included it here. The high P values 
of the three Galactic field MSPs listed above: B1937+21, 
J0218+4232 and PSR J1841+0130, imply that these MSPs 
are quite young. We discuss their true ages in Sections 18.3.31 
and [8X41 



4-3.1 Formation of a submillisecond pulsar 

The possible existence of submillisecond pulsars (P < 1 ms) 
is a long standing and intriguing question in the literature 
since it is importa nt for constraining the eq uation-of-state of 
neutron stars (e.g. lHaensel fc Zdunik 19891). There has been 
intense, but so far unsuccessful , efforts to detect these ob- 
jects in either radio or X-rays |D'Amicdl2000l : iKeith et al.l 
l20ld ; |Patrunoll20ld ). We note from equation © that if the 
inner edge of the accretion disk (roughly given by r mag ) ex- 
tends all the way down to the surface of the neutron star we 
obtain the smallest possible value for P cq ~ 0.6 ms (depend- 
ing on M and assuming ui c = 1). It is possible to achieve 
r mag ~ R from a combination of a large value of M and a 
small value of B, cf. equation ((2j. 



4.4 Amount of mass needed to spin up a pulsar 

4-4-1 Accretion-induced magnetic field decay 

The widely accepted idea of accretion-induced magnetic 
field decay in neutron stars is based on observational ev- 
idence (e.g. iTaam fe van den Heuvellll9 86: Shibazak i et al.l 
1 19891 ; IVan Den Heuvel fc Bitzarakil 1 1994! ) During the recy- 
cling process the surface B-field of pulsars is reduced by 
several orders of magnitude, from values of 10 11 " 12 G to 
10 7-9 G. However, it is still not understood if this is caused 
by spin- down induced flux exp ulsion of the core proton 
fluxoids (|Srinivasan et all Il990l ). or if the B-field is con- 
fined to the crustal regions and decays due to diffusion 
and Ohmic dissipation, as a result of a decreased electri- 
cal conductivity when heating eff ects set in from nuclear 
burning of the accreted ma terial (|Geppert fc Urpml 1 19941 : 
iKonar fc Bhattacharva|[T9 97). or if the field decay is simply 
caused by a di amagnetic screening by the accreted plasma 
- see review bv lBhattacharval l|2002i ) and references therein. 
Although there has been attempts to model or empirically fit 
the magnetic field evo l ution of accreting neutron stars (e.g. 
Shibazaki etafl 1 19891 : IZhang fc Koiimai I2OO6I : IWang et all 
20111) ihe results are quite uncertain. One reason is that it is 
difficult to estimate how much mass a given recycled pulsar 
has accreted. Furthermore, we cannot rule out the possi- 
bility that some of these neutron stars may originate from 
the accretion-induced collapse of a massive white dwarf, in 
which case they might be formed with a high B-field near 
the end of the mass-transfer phase. 

To model the B-field evolution for our purpose, which is 
to relate the spin period of pulsars to the amount of mass 
accreted, we make the following assumptions: 

• The B-field decays rapidly in the early phases of the 
accretion process via some unspecified process (see above). 

• Accreting pulsars accumulate the majority of mass 
while spinning at/near equilibrium. 

• The magnetospheric boundary, r mag is approximately 
kept at a fixed location for the majority of the spin-up phase, 
until the mass transfer ceases during the RLDP, cf. Sec- 
tion [S] 

• During the RLDP the B-field of an AXMSP can be 
considered to be constant since very little envelope material 
(~ 0.01 Mq) remains to be transfered from its donor at this 
stage. 
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4-4-2 Accreted mass vs final spin period relation 

The amount of spin angular momentum added to an accret- 
ing pulsar is given by: 



AJ* 



n(uj, t) M(t) yfGM{t)r mm (t)i(fil dt 



(11) 



where n{u) is a dimensionless torque, see Section [5] for a 
discussion. Assuming n(uj) = 1, and M(i), r mag (t) and £(£) 
to be roughly constant during the major part of the spin-up 
phase we can rewrite the expression and obtain a simple for- 
mula (see also equation [T]) for the amount of matter needed 
to spin up the pulsar: 



AM : 



2nl 



(12) 



Note, that the initial spin angular momentum of the pulsar 
prior to accretion is negligible given that Slo -C fl oq . To 
include all numerical scaling factors properly we can insert 
equations (fj?]l, © and © into equation (|12p . recalling that 
w = 4> ■ i~a, and we find: 




Equilibrium spin period, P gq (ms) 

Figure 6. The amount of mass needed to spin up a pulsar as a 
function of its equilibrium spin period using equation l|14jl . The 
different green curves correspond to various neutron star masses 
in steps of 0.2Mq, increasing upwards. The dashed red curve is 
the expression from equation l|15| l - sec text. 



AM, 



- _ j , --eg 

° q \G 2 M 



1/3 



(13) 



where f(a, £, <f>, lj c ) is some dimensionless number of order 
unity. Once again we can apply the relation between moment 
of inertia and mass of the neutron star (see Section I4.2.2p 
and we obtain a simple convenient expression to relate the 
amount of mass to be accreted in order to spin up a pulsar 
to a given (equilibrium) rotational period: 



AM eq = 0.22 M Q 



(M/Mq 



\l/3 



yt/3 



(14) 



assuming that the numerical factor f(a,£,(j),uu c ) = 1. 
In the above derivation we have neglected minor effects re- 
lated to release of gravitational binding energy of the ac- 
creted material - see e.g. e quation (22) i n lTauris fc Savoniie 
1999), and in particular iBeiger et al.l |201ll ') and Bagchi 



20111 ) for a more detailed, general discussion including var- 
ious equations of state, general relativity and the critical 
mass shedding spin limit. However, since the exchange of 
angular momentum takes place near the magnetospheric 
boundary the expression in equation (|14[) refers to the bary- 
onic mass accreted from the donor star. To calculate the in- 
crease in (gravitational) mass of the pulsar one must apply 
a reducing correction factor of ~ .85 — 0.90, depending o n 
the neutron star equation-of-state (|Lattimer fc Y ahil 1989). 

In Fig. [6]we show the amount of mass, AM eq needed to 
spin up a pulsar to a given spin period. The value of AM oq 
is a strongly decreasing function of the pulsar spin period, 
P e q. For example, considering a pulsar with a final mass of 
1.4 Mq and a recycled spin period of either 2 ms, 5 ms, 10 ms 
or 50 ms requires accretion of 0.10 Mq, 0.03 Mq, 0.01 Mq or 
1O~ 3 M0, respectively. Therefore, it is no surprise that ob- 
served recycled pulsars with massive companions (CO WD, 
O-Ne-Mg WD or NS) in general are much more slow rota- 
tors - compared to BMSPs with He WD companions - since 
the progenitor of their massive companions evolved on a rel- 
atively short timescale, only allowing for very little mass to 
be accreted by the pulsar, cf. Section [6] 



4-4-3 Comparison with other work 

The simple expres sion in equation (|14p was also found by 
lAI par et al.l (j 19821 ). Alternatively, one can integrate equa- 
tion (|12l) directly which yields the maximum spin rate (min- 
imum period) th at can be attained by a ne utron star accel- 
erated from rest l)Lipunov fc Postnovlll984l ): 



3nl 



(m 3/2 



M- 



3/2 



(15) 



AM cq is the mass of the neutron star 



where M ini t = M 
prior to accretion. In this expression it is still assumed that 
I and r mag remain constant during the accretion phase. For 
comparison, the expression above is also plotted in Fig.[6]as a 
dashed line assuming AM eq <C M, £ = 1 and r mag = 22 km 
(for example, if B = 10 s G, M = 0.1M Ed d,M = 1.4M , 
0=1, see equation ([2])). Although the latter expression is 
simplified with a chosen numerical value of r mag , and the two 
expressions differ by more than a factor of two if P > 10 ms, 
the overall match is fairly good. 

One the one hand, the value of AM oq should be re- 
garded as a lower limit to the actual amount of mate- 
rial required to be transfered to the neutron star, even 
at sub-Eddington rates, since a non-negligible amount 
may be ejected from the pulsar magnetosphere due 
to magnetodipole wave pr essure or the propeller effect 
(111 arionov fc Sunvaevl Il975f). Furthermore, a ccretion disk 
instabilities (|Pringld Il98ll : Ivan Paradiisl ll99fJ l are also re- 
sponsible for ejecting part of the transfered material . Re- 
cently it was demonstrated by I Antoniadis et al.l ( 2012) that 
the accretion efficiency in some cases is less than 50 per cent, 
even in short orbital period binaries accreting at sub- 
Eddington levels. On the other hand, we did not take into 
account the possibility of a more efficient angular momen- 
tum transfer early in the accretion phase where the value of 
r mag (the lever arm of the torque) could have been larger if 
the B-field did not d ecay rapidly. In thi s respect it is some- 
what surprising that IWang et alj (|201ll ) , using a model for 
accretion induced B-field decay, find that it requires more 
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mass to be accreted to obtain a given spin period (see their 
fig. 6) compared to our work. This could perhaps (partly) be 
related to the fact that Q. is small during the initial accretion 
phases. 



4.5 Spin-relaxation timescale 

Above we have calculated the amount of mass needed to spin 
up a pulsar to a given equilibrium spin period. However, we 
must be sure that the accretion-torque can actually trans- 
mit this acceleration on a timescale shorter than the mass- 
transfer timescale. To estimate the spin-relaxation timescale 
(the time needed to spin up a slowly-rotating neutron star 
to spin equilibrium) one can simply consider: ttorque = J/N 
where J = 2nl /P eq and TV = Msf GMr mag £ which yields: 

1/7 



4G 2 M 2 
„B 8 R 24 M 3 

50Myr B 8 " 8/ 



¥1 

M 



0.1 M E 



-3/7 



M 



1AM, 



(16) 

17/7 



or equivalently, ttorquc = AM/M, see equation (|12J1 . If 
the duration of the mass-transfer phase, tx is shorter than 
ttorquc (this is the case if the mass-transfer rate \M2\ S> M) 
then the pulsar will not be fully recycled. In the following 



section we examine the importance of t 



torque 



relative to the 



RLO turn-off timescale at the end of the mass-transfer pro- 
cess and how their ratio may affect the pulsar spin period. 



5 ROCHE-LOBE DECOUPLING PHASE - 
RLDP 

An important effect on the final spin evolution of accret- 
ing pulsars, related to the Roche-lobe decoupling phase 
(RLDP), has recently been demonstrated (|Taurisll2012r ). The 
rapidly decreasing mass-transfer rate during the RLDP re- 
sults in an expanding magnetosphere which may cause a 
significant braking torque to act on the spinning pulsar. 
It was shown that this effect can explain why the recy- 
cled radio MSPs (observed after the RLDP) are significantly 
slower rotators compared to the rapidly spinning accret- 
ing X-ray MSPs (observed before the RLDP) This differ- 
ence in spin periods was first noted bv lHessela |2008T ). How- 
ever, only for MSPs with B > 10 8 G can the difference in 
spin periods be partly understood from regular magneto- 
dipole and plasma current spin-down of the recycled radio 
MSPs. In this section we investigate the RLDP effect on 
IMXBs and LMXBs and demonstrate different outcomes 
for BMSPs with CO and He WD companions. The pur- 
pose of the computations is to follow the spin evolution 
of the accreting X-ray millisecond pulsar (AXMSP) when 
the donor star decouples from its Roche-lobe, and to cal- 
culate the initial spin period of the recycled radio pulsar 
once the mass-transfer ceases. We model the main effect (the 
RLDP which causes the magnetosphere to expand dramat- 
ically) on the general spin evolution and ignore other, less 
known or somewhat hypothetical, dynamical effects during 
this epoch, such as warped disks, disk-magnetosphere in- 
stabilities, variations in the mass transfer rate caused by 
X-ray irradiation effects of the donor star, or transitions be- 
tween a Keplerian thin disk and a sub-Keplerian, advection- 



dominated accretion flow (ADAF) (e.g. ISpruit fc Taam 


Il993|; [Nelson et al.l 19971: Yi et al. 


1997: van Kerkwiik et al. 


1998;: iLi & Wickramasinshe 199S 


: Locsei & Melatosl 12004!: 


iDai & Lill2006l; ICamero-Arranz et al.l 201 ll. and references 



therein). Although many of these effects may be less im- 
portant for the RLDP they could perhaps explain some 
of the frequent torq ue reversals observed in X-ray pulsars 
l|Bildsten et al.lll997h . 



5.1 Onset of a propeller phase 

In the final stages of the X-ray phase, when the donor star is 
just about to detach from its Roche-lobe, the mass-transfer 
rate decreases significantly. This causes the ram pressure of 
the incoming flow to decrease whereby the magnetospheric 
boundary, and thus the coupling radius r mag , moves out- 
ward. When this boundary moves further out than the coro- 
tation radius (r mag > r co ) given by: 



1/3 



(GM_\ 

\m ) 

17 km P 2/ s 3 



(17) 



M 



1.4 M, 



1/3 



a centrifugal barrier arises since the plasma flowing to- 
wards the neutron star couples to the field lines in a 
super-Keplerian orbit. The material is therefore presum- 
ably ej ected away from the neutr on star in this propeller 
phase (|lllarionov fc Sunvaevl ll975). This ejection of mate- 
rial causes exchange of angular momentum between the now 
relatively fast spinning neutron star and the slower spinning 
material at the edge of the inner disk. The result is a brak- 
ing torque which acts to slow down the spin of the pulsar. 
This causes the corotation radius to move outwards too and 
the subsequent evolution is determined by the rate at which 
these two radii expand relative to one another. However, 
during the RLDP M decreases too rapidly for r co to keep 
up with the rapidly expanding r mag and the equilibrium spin 
phase comes to an end. 

In our model we assumed, in effect, that the inner edge 
of the accretion disk, rdisk follows r mag and that the cen- 
trifugally expelled material has sufficient kinetic energy to 
be gravitationally unbound during the propeller phase when 
^mag > r co . These assumptions may not be en tirely correct 
i|Spruit fc Taaml 1 19931: iRappaport et al.l 120041 ) . For exam- 
ple, iD'Angelo fc Spruit! (|20ld . 1201 ll . I2012T) recently demon- 
strated that accretion disks may be trapped near the coro- 
tation radius, even at low accretion rates. In this state the 
accretion can either continue to be steady or become cyclic. 
Furthermore, based on energy considerations they point out 
that most of the gas will not be ejected if the Keplerian 
corotation speed is less than the local escape speed. This 
will be the case if rdisk < 1.26 r co . However, this trapped 
state only occurs in a narrow region around the corotation 
radius (Ar <g r co ) and given that r mag /r- co often reaches a 
factor of 3 — 5 in our models at the end of the RLDP we find 
it convincing that the RLDP effect of braking the spin rates 
of the recycled pulsars is indeed present (see iTaurisl (|2012T ) 
for further details). 
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5.1.1 The magnetosphere-disk interaction 

In our numerical calculations of the propeller phase we in- 
cluded the effect of additional spin-down torques, acting on 
the neutron st ar, due to both magne tic field drag on the 
accretion disk jRappaport et aLll2004T l as well as magnetic 
dipole radiation (see equation (|3|), although these effects 
are usually not dominant. It should be noted that the mag- 
netic stress in t he disk is also related to the critical f astness 
param eter, u c (jGhosh fe Lambl Il979bh . We follow iTaurij 
(|2012T ) and write the total spin torque as: 

iVtotai = n(w) (AVGMr mag £ + ^—^j - (18) 
where 

n(w) = tanh \JJ^-J ( 19 ) 

is a dimensionless function, depending on the fastness pa- 
rameter, uj = n^/nicO-mag) = (^mag/rco) 3/2 , which is in- 
troduced to model a gradual torque change in a transition 
zone near the magnetospheric boun dary. The width of thi s 
zone has been shown to be small l|Spruit fc Taaml 1 1993ft . 
corresponding to i5 u < 1 and a step-function-like behav- 
ior n(u>) = ±1. In our calculations presented here we used 
Su = 0.002, f = 1, rdi s k = fmag and also assumed the mo- 
ment of inertia, I, to be constant during the RLDP. The 
latter is a good approximation since very little material is 
accreted during this terminal stage of the RLO. 

5.1.2 Radio ejection phase 

As the mass-transfer rate continues to decrease, the mag- 
netospheric boundary eventually crosses the light-cylinder 
radius, ri c of the neutron star given by: 

Tic = c/O, ~ 48 km P mB (20) 

When r mag > ric the plasma wind of the pulsar can stream 
out along the open field lines, providing the necessar y condi- 
tion for the radio emiss ion mechanism to turn on fe.g. lMichell 
1 199 ll : ISpitkovskvll2008l ) . Once the radio millisecond pulsar is 
activated the presence of the plasma wind may prevent fur- 
ther accretion from th e now weak flow of material at l ow M 
jKluzniak et al.ll 19881 ). As discussed bv lBurderi et all |200ll . 
2002) this will be the case when the total spin-down pressure 
of the pulsar (from magneto-dipole radiation and the plasma 
wind) exceeds the inward pressure of the material from the 
donor star, i.e. if: _E ro t/(47rr 2 c) > Pdisk- However, it is inter- 
esting to notice that t he recycled pulsar PSR J1023+0038 
( Archibald ct al. 2009]) recently (within a decade) turned on 
its radio emission and it would indeed be quite a coincidence 
if we have been lucky enough to catch that moment - unless 
recycled pulsars evolve through a final phase with recurrent 
changes between accretion and radio emission modes. 

5.2 Slow RLDP 

In Fig. [7] we have shown a model calculation of the RLDP of 
an LMXB with an original donor star mass of 1.1 Mq. The 
outcome was the formation of a BMSP with P — 5.2 ms (as- 
suming B — 1.0 x 10 8 G) and a He WD companion of mass 



0.24 Mq , orbiting with a period o f 5.0 days. It has b e en ar - 
gued by ITauris fc Savoniid \ 19991 ) : lAntoniadis et ail (|2012T l 
that the majority of the transfered material in some LMXBs 
is lost from the system, even for accretion at sub-Eddington 
rates. Therefore, we assumed an effective accretion efficiency 
of 30 per cent in our LMXB model, i.e. M = 0.30 |M 2 |, where 
IM2I is the RLO mass-transfer rate from the donor star. The 
accreted^ mass-transfer rate M is shown in the upper left 
panel. 

The three phases: equilibrium spin (corresponding to r mag ~ 
r co ), propeller phase (r co < r mag < ric) and radio pulsar 
(fmag > 7ic) ar e clearly identified in this figure. During the 
equilibrium spin phase the rapidly alternating sign changes 
of the torque (partly unresolved on the graph in the lower 
left panel) reflect small oscillations around the semi-stable 
equilibrium, corresponding to successive small episodes of 
spin up and spin down. The reason is that the relative lo- 
cation of r mag and r co depends on the small fluctuations in 
M. Despite applying an implicit coupling scheme in the code 
and ensuring that our time steps during the RLDP (about 
8 x 10 4 yr) are much smaller than the duration of the RLDP 
(~ 10 8 yr), our calculated mass-transfer rates are subject 
to minor numerical oscillations. However, these oscillations 
could in principle be physical within the frame of our simple 
model. Examples of physical perturbations that could cause 
real fluctuations in M, but on a much shorter timescale, 
are accretion disk instabilities and clumps in the transfered 
material from the donor star (our donor stars have active, 
convective envelopes during the RLDP). 

At some point the equilibrium spin is broken. Initially, 
the spin can remain in equilibrium by adapting to the de- 
creasing value of M. However, further into the RLDP the 
result is that r mag increases on a timescale faster than the 
spin-relaxation timescale, ttorquc at which the torque can 
transmit the effect of de celeration to the neutron star, and 
therefore r mag > r co fsee lTauris]|2012l . and supporting online 
material therein). During the propeller phase the resultant 
accretion torque acting on the neutron star is negative, i.e. 
a spin-down torque. When the radio pulsar is activated and 
the accretion has come to an end, the spin-down torque is 
simply caused by the magneto-dipole radiation combined 
with the pulsar wind (see Section 223} • 

The aim here is to calculate the initial spin period of the 
recycled radio pulsar, Po once the mass-transfer ceases. In 
this LMXB model calculation the value of the spin period in- 
creased significantly during the RLDP. The reason for this is 
that the duration of the propeller phase (i.e. the RLDP) is a 
substantial fraction of the spin-relaxation timescale, ttorquc- 
Using equation (|16p we find ttorquc = 195 Myr whereas the 
RLDP lasts for £r.ldp — 56 Myr which is a significant frac- 
tion of ttorquc (tRLDp/ttorquc = 0.29). Therefore, the RLDP 
has quite a significant impact on the spin evolution of the 
neutron star. It is seen that the spin period increases from 
3.7 ms to 5.2 ms during this short time interval, i.e. the 
pulsar loses 50 per c ent of i t s rota tional energy during the 
RLDP. As shown by iTaurisl (|2012T ) this RLDP effect is im- 
portant for understanding the apparent difference in spin 



1 Strictly speaking this is the estimated mass-transfer rate re- 
ceived from the inner edge of the accretion disk. Some of this 
material is ejected during the propeller phase. 
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Age (Gyr) Spin period of pulsar, P (sec) 

Figure 7. Example of a slow Roche-lobe decoupling phase (RLDP) and the final spin evolution of a pulsar formed with a He WD 
companion in an LMXB. The left side panels show: The accretion rate of the pulsar, M as a function of the total age of the donor star 
(top); the evolution of the location of the magnetospheric boundary, r mag , the corotation radius, r co and the light-cylinder radius, r\ c 
(centre); and the resulting braking torque (bottom). The right side panels show: The evolution of the pulsar spin period as a function of 
the total age of the donor star (top); the evolution of the pseudo characteristic age of the pulsar (centre); and the RLDP evolutionary 
track leading to the birth location in the PP-diagram shown as a blue diamond (bottom). The three epochs of the RLDP: equilibrium 
spin, propeller phase and radio pulsar are easily identified in the left side panels. The propeller phase is grey shaded. In these model 
calculations we assumed sin a = <^> = £ = ai c = l and <5„ = 0.002. It is interesting to notice how the changes in the accretion rate, M affects 
the magnetospheric coupling radius, r mag which again affects the spin-down torque and the equilibrium spin period, cf. equations J2j|, 
UJ and J7J. The reason why the spin of the pulsar (blue line, upper right panel) decouples from its equilibrium spin (black line) is that 
the spin-down torque cannot transmit the deceleration on a timescale short enough for the pulsar to adapt to its new equilibrium. In 
this example, the duration of the RLDP is a significant fraction of the spin-relaxation timescale of the accreting neutron star and thus 
the RLDP has a momentous impact on the rotational evolution. The birth spin period, Pq of this recycled pulsar is seen to deviate 
substantially from its original P cq calculated at the onset of the RLDP. Furthermore, the effect on the characteristic age for this BMSP 
is also quite significant. During the 56 Myr of the RLDP it doubles to become r = 13 Gyr at birth. Hence, many MSPs which appear to 
be old according to their r may in fact be quite young. See text for further discussions. 
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Figure 8. Example of a rapid Roche-lobe decoupling phase (RLDP) in an IMXB using our mass-transfer progenitor model of 
PSR J1614— 2230. The spin evolution of the pulsar is shown during the RLDP when the IMXB donor star terminated its mass transfer. 
For a description of the various panels and their labels, see Fig. [7] In the PP— diagram in the lower right panel the birth location, 
immediately following the RLDP, of the recycled radio pulsar is shown as a blue diamond for sin a = </> = £ = aj c = l. The red star 
indicates the presently observed location of PSR J1614— 2230. In Section I9.ll we argue that PSR J1614— 2230 is more likely to have 
evolved with sine* < 1, (j> > 1 and u c < 1 which results in a birth location much closer to its present location in the PP-diagram. In 
IMXBs where the donor star evolves on a relatively short timescale, the RLDP is rapid. Therefore, the spin evolution of the accreting 



X-ray pulsar (see blue line in upper right panel) freezes up near the original value of P 0< 



the recycled radio pulsar is born with a 



spin period, Pq which is close to the equilibrium spin period of the accreting X-ray pulsar calculated at the time when the donor star 
begins to decouple from its Roche-lobe. In this IMXB model the mass-transfer rate, M > M^d an d the accretion onto the neutron star 
is therefore limited by the Eddington mass-accretion rate (red dot-dashed line in the upper left panel). This explains the why the spin 
of the pulsar, P = P eq ~ const, (top right panel) until M < M^dy as opposed to the situation in Fig. [7] 
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period distrib utions between AXMSPfl and radio MSPs 
(|Hesseldl2008h . The black line labelled "P cq " in the upper 
right panel reflects the evolution of P in case the neutron 
star was able to instantly re-adjust itself to the changing 
P eq (equation [7} during the RLDP. However, the neutron 
star possesses a large amount of rotational inertia and such 
rapid changes in P are not possible to transmit to the neu- 
tron star given the limited torque acting on it. Therefore the 
spin evolution terminates at the end of the blue line labelled 
"Po". The small fluctuations of P oq during the equilibrium 
spin phase reflect fluctuations in M(t). 

The RLDP effect also causes a large impact on the char- 
acteristic age, to, of the radio pulsar at birth. In the central 
right panel we show how to doubles in value from 6.7 Gyr to 
13 Gyr during the RLDP (see also the solid line ending at the 
blue diamond in the PP-diagram in the lower right panel). 
The characteristic age at birth is given by: to = Po/2Po, 
where Po can be estimated from the assumed B-field of the 
pulsar at this epoch. Actually, the "characteristic age" of 
the pulsar during the RLDP (which is plotted in the central 
right panel) is a pseudo age which is only relevant at the end 
point (to) when the accretion has stopped. In order to con- 
struct the track leading to to we assumed B 2 oc PP = const. 
from equation @, leading to t oc P 2 . During the RLDP it 
is a good approximation to assume a (final) constant value 
of the neutron star surface magnetic flux density, B since 
99 per cent of the accretion onto the neutron star occurs 
before the final termination stage of the RLDP. In this case 
we assumed a final B-field strength of B = 1 X 10 s G - a 
typical value for BMSPs with He WD companions. 
This example clearly reflects why characteristic ages of 
MSPs are untrustworthy as true age ind icators s ince h ere the 
MSP is born with t ~ 13 Gyr (see also lTaurisI l|2012D V This 
fact is important since the characteristic ages of recycled pul- 
sars are often compared to the cooling ages of their white 
dwar f companions (e.g. see discussion in lvan Kerkwiik et al.l 
l2005h . 

5.3 Rapid RLDP (P cq freeze- up) 

In Fig. [5] we have plotted the RLDP during the final stages 
of the IMXB phase AB using our Case A model calcula- 
tion of PSR J1614-2230 presented in Paper I. The initial 
conditions were a 4.5 Mq donor star in an IMXB with an 
orbital period of 2.20 days. The final system was a BMSP 
with a 0.50 M CO WD orbiting a 1.99 M recycled pulsar 
with an orbital period of 8.7 days - resembling the char- 
acteristics of the PSR J1614— 2230 system. The accretion 
rate onto the neutron star, M (see blue line in upper left 
panel) is Eddington limited (red dot-dashed line) and hence 

M — min ( \Mi\, Medd) where \M%\ represents the mass- 
transfer rate from the donor star (purple line). When the 
age of the donor star is ~ 170.0 Myr the mass-transfer rate 

2 The majority of the AXMSPs have very small P or b (typically 
1-2 hours) and substellar companions (< 0.08 Mq) and hence 
most of these AXMSPs are not true progenitors of the general 
population of radio MSPs. However, the important thing is that 
the pulsars are spun up to rotational periods which seem to be 
almost independent of orbital period up to 200 days, as seen in 
the central panel of Fig. [3] 



becomes sub- Eddington and M decreases with \M?\. From 
this point onwards r mag increases as expected, as a result of 
the decreasing ram pressure, and it increases at a rate where 
r co cannot keep up with it (see central panel, left column) 
and therefore the pulsar enters the propeller phase. After 
less than 0.3 Myr the accretion phase, and thus the pro- 
peller phase, is terminated when r mag > ri c and the radio 
ejection mechanism is activated. 

During the equilibrium spin phase the net accretion 
torque acting on the pulsar is close to zero. This is not only 
due to our transition zone approximation (equation (|19l) ~) 
where the torque vanishes when r mag ~ r co . If the transition 
zone was arbitrary thin (8^ -^g 1), and n(uj) would be a step- 
like function going from +1 to —1 (i.e. switch mode), there 
would be rapidly alternating sign changes of the torque dur- 
ing the equilibrium spin phase. Whereas small fluctuations 
in M would lead to such rapid torque oscillations for sub- 
Eddington accretion (as seen in Fig. [7} the sharply defined 
Eddington limit applied here prevents such fluctuations dur- 
ing the equilibrium spin phase. However, the time-averaged 
torque would be close to zero in any case. 

In this IMXB model calculation the value of the spin pe- 
riod only increased very little during the RLDP. The reason 
for this is that the duration of the RLDP (t RL DP ^ 0.28 Myr) 
is relatively short compared to the spin-relaxation timescale 
(ttorquc = 6.8 Myr). The small ratio tRLDp/itorquc — 0.04 at 
the onset of the RLDP causes the spin period to "freeze " 
at the original value of P eq (see also iRuderman et al.lll98"9T ). 
Note, the values quoted above as well as rough values of 
^torque = (2ir/P)I/N and r mag ~ r co ~ 22 km can be 
checked by reading numbers from the plots and using equa- 
tions (J2J) and (TTTJ). The relevant numbers are: B = 5 x 10 s G 
(a typical value for BMSPs with CO WD companions), 
M ~ M E dd, M ~ 1.99 M e , P oq = 1.28 ms and logN = 34.8. 
For the spin-up process discussed here we assumed again 
sin a = <f> = £ = (j a = l. 

5.3.1 Rapid RLDP with a helium star donor 

As discussed in Section [7] (and seen in Fig. I10[) helium 
donor stars evolve on a short nuclear timescale compared 
to a normal hydrogen-rich donor. This leads to a rapid 
accretion turn-off lasting only a few 10 4 yr and therefore 
tRLDp/ttorquo <C 1, which also in these systems causes Po to 
"freeze" at the original value of P cq . 

5.4 RLDP in IMXB vs LMXB systems 

The major difference between the RLDP in an IMXB (of- 
ten leading to a BMSP with a CO WD) and the RLDP 
in an LMXB (leading to a BMSP with a He WD) is the 
time duration of the RLDP relative to the spin-relaxation 
timescale, i.e. the tRLDp/ttorque-ratio. The mass-transfer in 
X-ray binaries proceeds much faster in IMXBs compared to 
LMXBs. Hence, also the termination of the RLDP is shorter 
in IMXBs, and therefore the irtLDp/itorquc-ratio is smaller, 
compared to LMXBs. The conclusion is that, in general, we 
expect a significant RLDP effect only in LMXBs. (It may be 
possible, however, that the RLDP effect could play a more 
significant role in some IMXB systems if they leave behind a 
neutron star with a relatively high B-field, thereby decreas- 
ing ttorquc)- 
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Spin period, P (sec) 

Figure 9. The observed spin period distribution of recycled pul- 
sars in the Galactic disk with P < 100 ms for pulsars with He WD 
and CO WD companions. The median spin periods for these pop- 
ulations are 3.9 ms and 16 ms, respectively. For comparison we 
have also shown the spin distribution of those recycled pulsars 
with neutron star companions and P < 100 ms. 



6 OBSERVED SPIN-PERIOD DISTRIBUTIONS 

We now investigate if the empirical pulsar spin data can be 
understood in view of the theoretical modelling presented in 
the last couple of sections. 



6.1 Recycled pulsars with WD companions 

In Fig. [5] we have plotted histograms of the observed spin 
period distribution of recycled pulsars with He WD and 
CO WD companions. There are 57 known recycled pulsars 
in the Galactic disk with a He WD companion and a spin pe- 
riod P < 100 ms (with or without a measured P) and these 
pulsars have a median spin period of 3.9 ms (the average spin 
period is 10 ms). The median spin period of the 16 known 
recycled pulsars with a CO WD companion and P < 100 ms, 
however, is 16 ms (the average value is 23 ms). All observed 
BMSPs have lost rotational energy due to magneto-dipole 
radiation and a pulsar wind since they appeared as recy- 
cled pulsars. Hence, their current spin period, P is larger 
than their P aq . Unfortunately, we have no empirical con- 
straints on the spin-down torque acting on recycled radio 
MSPs and therefore their braking index remains unknown, 
cf. Section [5] Furthermore, in Section [5] we demonstrated 
that BMSPs formed in LMXBs lose rotational energy dur- 
ing the RLDP. Correcting for both of these spin-down effects 
we here simply assume P = v2P eq . The resulting median 
values of P Gq are thus found to be 2.8 ms and 11 ms, re- 
spectively. For a 1.4 Mq neutron star these median equi- 
librium spin periods correspond to AM eq = 0.06 Mq and 
AM eq = 0.01 Mq for pulsars with a He WD or a CO WD 
companion, respectively. In other words, to explain the ob- 
served period distribution of these two classes of recycled 
pulsars we conclude that those pulsars with CO WD com- 
panions typically have accreted less mass by a factor of 6, 
compared to the pulsars with He WD companions. This con- 
clusion remains valid for other scalings between the present 
spin period, P and the initial spin period, P cq - even for 



assuming P — P cq . We conclude from Fig. [9] that BM- 
SPs with CO WDs usually have evolved via early Case B 
or Case C RLO (the latter leading to CE evolution). The 
reason is that these two paths involve mass transfer on a 
much shorter timescale, and hence lead to smaller amounts 
of transfered mass, leading to less effective recycling and 
thereby longer spin periods, as seen among the observed 
pulsars in Fig. [9] (see also Section [7J. Furthermore, the BM- 
SPs with CO WDs also have relatively large values of P. 
This reflects larger values of B, which is expected since less 
mass was accreted in these systems, cf. Section 14.4.11 The 
only exception known from the description outlined above 
is PSR J1614-2230 which is discussed in Section H 
For BMSPs with He WD companions, we notice from the 
orbital period distribution in Fig. [3] (central panel) that the 
majority of these systems seem to have evolved through 
Case B RLO with an initial LMXB orbital period of P OT b > 
ldav ( | Tauris fc Savoniiei 1 19991 : IPodsiadlowski et~al ] I2OO2I : 
lTaurisll201ll . and references therein). 



6.2 Double neutron star systems 

In Fig. [9] we have also plotted the spin distribution of 
recycled pulsars with neutron star companions and P < 
100 ms. These syste ms originate from high-ma s s X-ra y bina- 
ries (HMXBs), e.g. iTauris fc van den Heuvell (|2006T l. Their 
spin distribution is skewed to larger periods compared to 
the recycled pulsars with CO WD companions. This is ex- 
pected since their massive stellar progenitors evolved even 
more rapidly than the donor stars of the IMXBs. It is gen- 
erally believed that double neutron star systems evolved via 
a CE and spiral-in phase since the mass tran sfer in HMXB s 
is always dynamic ally unstable - see however iBrownl (|l995h ; 
iDewi et al. (|2006h for an alternative formation mechanism 



via a double-core scenario. 



6.3 Formation of recycled pulsars with 

Porb > 200 days 

As mentioned earlier (see also Fig. [3}, pulsars with P or b > 
200 days all have slow rotational spins. This fact could be 
explained if they had accreted little mass. If their progeni- 
tors were wide-orbit LMXBs they had experienced continu- 
ous super-Eddingto n mass transfer on a sh ort timescale of 
the order 10 Myr l|Tauris fc Savoniid[l999l ). However, even 
if their accretion rates had been restricted to the Edding- 
ton limit these accreting pulsars could still have received up 
to 0.3 M©. This amount is 10 times larger than needed to 
spin up a pulsar to 5 ms according to equation (|14|l . There- 
fore, to explain the slow spin periods of these wide-orbit 
recycled pulsars one may speculate that the accretion pro- 
cess was highly inefficient due to, for example, enhanced 
accretion disk instabilities (Ivan Paradiisllili 6; D ubus et al.l 
1999) . This could be related to the advanced evolutionary 
stage of their donor stars on the red giant branch (RGB). 
These stars had deep convective envelopes which may lead 
to enhanced clump formation in the transfered material. 
Alternatively, one may ask if these systems descended from 
wide-orbit IMXBs which underwent a CE and a mild spiral- 
in. In this case their companion stars should be fairly mas- 
sive CO WDs (^ 0.5 Mq, see fig. 1 in Paper I). However, 
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quite a few of these wide-orbit recycled pulsars seem to 
have relatively low-mass He WDs accordin g to the i r mas s 
functions. In fact it wa s pointed ou t by Taurisl i|l996l ); 
iTauris fc Savoniid (|l999h ; IStairs et al] (|2005r ) that some of 
those WD masses may even be significantly less massive 
than expected from the (Mwd, P rt>)-relation which applies 
to He WDs descending from LMXBs, although this analysis 
still has its basis in small numbers. It is therefore impor- 
tant to observationally determine the masses of the WDs in 
these wide-orbit pulsars systems and settle this evolutionary 
question. 



7 SPINNING UP POST COMMON ENVELOPE 
PULSARS 

An interesting question is if our derived relation between 
pulsar spin and accreted mass (equation I14[) can also ex- 
plain the observed spin periods of the BMSPs with mas- 
sive WD companions and short orbital periods, i.e. sys- 
tems we expect to have evolved through a CE evolution. 
Let us consider th e important system P SR J1802-2124 re- 
cently observed bv lFerdman et all (2010). The nature of its 
0.78 ± 0.04 Mq CO WD companion, combined with an or- 
bital period of 16.8 hr and a post-accretion pulsar mass of 
only 1.24 ± 0.11 Mq makes this one of our best cases for 
a system which evolved via a CE and spiral-in phase. The 
present spin period of the pulsar is 12.6 ms. Assuming an 
initial equilibrium spin period of 10 ms we find from equa- 
tion (|14[) that AAf cq ~ 0.01 Mq which is a much more ac- 
ceptable value for a system undergoing CE-evolution (com- 
pared to the 0.10 Mq needed for spinning up a BMSP with 
a He WD companion to a spin period of a few ms). However, 
one still has to account for the 0.01 Mq accreted. It is quite 
likely that this accretion occurred after the CE-phase and 
we now consider this possibility. 



7.1 Case BB RLO in post common envelope 
binaries 

The total duration of the CE-phase itself is probably so 
short (< 10 3 yr, e.g . |Podsiadlowsklll200ll ; iPassv et al]|2012t 
llvanova et al. 1 120121 ) that no substantial accretion onto the 
neutron star is possible with respect to changing its rota- 
tional dynamics - even if the envelope of the donor star is 
not ejected from the system until long after the in-spiral of 
the secondary star which occurs on a dynamical timescale 
of a few P or b- Thus we are left with two possibilities: 1) the 
recycling of the pulsar occurs during the subsequent stable 
RLO of a naked helium star following the CE (Case BB 
RLO), or 2) the recycling occurs due to wind accretion from 
the naked helium or carbon-oxygen core left behind by the 
donor star after it loses its envelope in the CE-phase. 
In Fig. [TO] we show one of our model calculations of Case BB 
RLO. In our example we assumed a 1.40 A/q helium star 
(Y — 0.98) as the donor and a 1.35 Mq neutron star as the 
accretor. The giant phase of a helium star is short lived and 
in this case the mass transfer lasts for about 170.000 yr. 
(This time interval, however, is much longer than the CE- 
event.) Assuming the accretion onto the neutron star to be 
Eddington limited will therefore lead to an accretion of only 
about 7x 10 -3 Mq. Although this is a small amount, it is still 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
1.40 M Q He — star donor 




1.35 M Q neutron star 




P init = °* 50 da V s 




AM = 0.002 M Q 






M Edd 







Age (Myr) 

Figure 10. A model calculation of Case BB RLO. The plot 
shows mass-transfer rate as a function of age for a 1.4 Mq helium 
star donor. The initial orbital period is 0.50 days and the neutron 
star mass is 1.35 Mq. After the RLO the donor star settles as a 
0.89 Mq CO WD orbiting the recycled pulsar with an orbital 
period of 0.95 days. The mass-transfer phase only lasts for about 
170.000 yr. Even if the neutron star only accretes ~ 0.002 Mq 
(assuming an accretion efficiency of just 30%) it is still able to 
spin up the pulsar to a spin period of 36 ms. If the neutron star 
accretes all matter at the Eddington limit it can be spun up to 
14 ms. A comparison can be made with the observed spin periods 
of BMSPs with CO WD companions in the bottom panel of Fig. [3] 



sufficient to spin up the recycled pulsar to a period of 14 ms, 
which is even faster than the observed median spin period 
of BMSPs with CO WD companions (see Fig. [9j). Assum- 
ing the accretion efficiency to be only 30 per cent of M^dd 
would still spin up the pulsar to 36 ms. A helium donor star 
with a lower mass can even spin up the pulsar below 14 ms. 
We found that for a 1.10 Mq helium star donor the neutron 
star can accrete 50 per cent more (resulting in P cq = 11ms) 
due to the increased nuclear timescale of the shell burn- 
ing of a lower mass helium star. Hence, based on rotational 
dynamics we conclude that Case BB RLO is a viable forma- 
tion channel to explain most of the observed BMSPs with 
CO WD companions. A more systematic study of Case BB 
RLO for various system configurations is needed in order 
to find the minimum possible spin period of BMSPs with 
CO WD companions formed via this formation channel. 

It should be noted that a pulsar may also be recycled 
via Case BA RLO, i.e. from RLO which is initiated while 
the companion st ar is still on the helium main sequence (see 
iDewi et al.ll2002l . for detailed calculations). In this case the 
mass-transfer phase will last up to ten times longer, com- 
pared to Case BB RLO, providing more than sufficient ma- 
terial to recycle the pulsar effectively. However, in order to 
initiate RLO while the helium star donor is still burning core 
helium on its main sequence requires a tight, fine-tuned in- 
terval of orbital separation following the common envelope. 
Typical orbital periods needed are 1-2 hours. Since Case BB 
RLO is initiated for a much wider interval of larger orbital 
periods (up to several tens of days) we expect many more 
systems to evolve via Case BB RLO compared to Case BA 
RLO. 
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7.2 Wind accretion prior to Case BB RLO 

Finally, we tested if the neutron star could acquire any 
significant spin up from wind accretion in the epoch be- 
tween post-CE evolution and Case BB RLO. To produce 
CO or ONeMg WD remnants via Case BB RLO one 
would expect a typical helium star mass of 1.1 — 2.2 Mq 
prior to the mass transf er. These star s have luminosities of 
log(L/I/0) = 2.5 — 3.3 (|Langeilll989] ) and typically spend 
less than 10 Myr on the He-ZAMS before they expand and 
initiate Case BB RLO. Their ma ss-loss rates are of the or - 
der M wind ~ lO -10 * ' 6 M© yr" 1 jjefferv fc Hamannll2010l ). 
Therefore, even if as much as 10 per cent of this ejected wind 
mass was accreted onto the neutron star it would accrete a 
total of at most ~ 10 -4 Mq prior to the RLO. We there- 
fore conclude that wind accretion prior to Case BB RLO is 
negligible for the recycling process of BMSPs with CO WD 
companions. 



7.3 Will the progenitor binaries make it to 
Case BB RLO? 

Before drawing conclusions about the Case BB RLO we 
should verify that the progenitor binary survives to follow 
this path. In Paper I we argued Case BB RLO does not work 
for PSR J1614-2230 because the envelope of the CO WD 
progenitor star is too tightly bound to be ejected during the 
spiral-in phasqfl I s that the case for all BMSPs? The an- 
swer is no. We find that if the post-CE orbital period (i.e. 
the observed P or b of a BMSP) is about 0.5 days or less then 
sufficient orbital energy will be released during spiral in to 
energetically allow for a successful ejection of the envelope 
of the original 5 — 7 Mq red giant donor star. The total 
(absolute) binding energy of the envelope of such a star is 
1.3 — 2.2 x 10 48 ergs which is equivalent to an orbital energy 
corresponding roughly to P or b ^ 0.5 days. Since we must re- 
quire Ai5 orb > -Benv to avoid a merger we therefore conclude 
that Case BB RLO is in principle a viable formation chan- 
nel for BMSPs with CO WD companions if these systems 
have short orbital periods. As argued earlier we consider 
PSR J1802-2124 (which has P orb = 16.8 hours) as a strong 
case for a BMSP which evolved though a common envelope. 
Whether or not one can finetune the Case BB RLO evo- 
lution to account for the recycling of this pulsar, or if the 
recycling was due to wind accretion from the exposed post- 
CE CO core of an AGB star or, thirdly, if accretion-induced 
collapse of an ONeMg WD was at work, is an interesting 
question to address. 



8 TRUE AGES OF MILLISECOND PULSARS 

Knowledge of the true ages of recycled radio pulsars is im- 
portant for comparing the observed population with the 
properties expected from the spin-up theory outlined in Sec- 
tion [4] All radio pulsars lose rotati onal energy with age 
and t he braking index, n is given by ^Manchester fe TavloJ 
1 19771 ): 



3 We even considered a conservative Case BB RLO evolution 
which is far from valid given that the mass-transfer rate is highly 
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Figure 11. Isochrones of pulsar evolution for a pulsar with known 
initial (Po, Pq) for different values of a constant braking index, n. 
Along any given isochrone n increases continuously from n = 1 
(top) to n = 15 (bottom). The initial position of the recycled 
test pulsar is shown by a red star and the dashed lines show 
evolutionary tracks for n = 2, 3 or 5, respectively. The dotted line 
yields a characteristic age, r = P/ (2P) of 12 Gyr. This line would 
intersect with the 12 Gyr isochrone exactly at n = 3 if Po <S P, 
which is not completely fulfilled. The dot-dashed isochrone line 
is a general solution to equation (1 23 P for n = 3, t = 6 Gyr and 
Pq = 2.5 ms, but unknown Pq. 



(21) 

Q.Q,/tl 2 . This deceleration 



Cl oc — fi n 

which yields (for n constant): 
law can also be expressed as: P oc P 2 ~ n and hence the slope 
of a pulsar evolutionary track in the PP-diagram is simply 
given by: 2 — n. Depending on the physical conditions under 
which the pulsar spins down n can take different values as 
mentioned earlier. For example: 



gravitational wave radiation n — 5 

B-decay or alignment or multipoles n > 3 

perfect magnetic dipole n = 3 

B-growth/distortion or counter-alignment n < 3 



(22) 



The combined 
rent spin-down 



and plasma cur- 
result in n 7- 3 
A simple integra- 
tion of equation (|21|l for a constant braking index (n 7- 1) 
yields the well-known expression: 



magnetic dipole 
torque m a y als o 
l|Contopoulos fc Spitkovskvl [2006). 



(n-l)P 



1- ^ 



(23) 



super-Eddington, see Fig. 1 101 above. 



where t is the so-called true age of a pulsar, which had an 
initial spin period Pq at time t = 0. 



8.1 Isochrones in the PP— diagram 

Unfortunately, one cannot directly use the equation above 
to obtain evolutionary tracks in the PP-diagram (even un- 
der the assumption of a constant n and a known initial 
spin period, Po). For chosen values of t, n and Po one can 
find a whole family of solutions (P, P) to be plotted as an 
isochrone, see the dot-dashed line in Fig. [11] (and see also 
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Figure 12. Evolutionary tracks of characteristic ages, t calculated as a function of true ages, t for recycled pulsars with a constant 
braking index of n = 2, n = 3, and n = 5. In all cases for n we assumed a constant M = IAMq, since = 1, a constant moment of inertia, 
I and for n = 5 we also assumed a constant cllipticity, e ^ 0. In the left panel we assumed in all cases an initial spin period of Po = 3.0 ms 
and varied the value of the initial surface magnetic flux density, Bo- In the right panel we assumed in all cases Bo = 10 8 G and varied 
Po- The dotted line shows a graph for r = t and thus only pulsars located on (near) this line have characteristic ages as reliable age 
indicators. We notice that recycled pulsars with small values of Po, resulting from either small values of Bo (left panel) and/or large 
values of Po (right panel), tend to have r>t, even at times exceeding the age of the Universe. In all degeneracy splittings of the curves 
the upper curve (blue) corresponds to n = 5, the central curve (red) corresponds to n = 3 and the lower curve (green) corresponds to 
n = 2. The asymptotic behavior of the curves can easily be understood from equation {25} — see text for explanations. 



iKiziltan fc ThorsettJlioioh . However, there is only one point 
which is a valid solution for a given pulsar and the above 
equation does reveal which point is correct. The problem is 
that the variables in equation (|23p are not independent and 
we do not know a priori the initial spin period derivative, Po. 
The evolution of the spin period is a function of both t, n, 
the initial period and its time derivative, i.e. P(t, n, Po, Po). 
Here, for simplicity, we assume n to be constant. To de- 
termine Po one must make assumptions about the initial 
surface magnetic flux density, Bo, the initial magnetic incli- 
nation angle, ao and the initial ellipticity, £o (the geomet- 
ric distortion which is relevant for rotational deceleration 
due to gravitational wave radiation). Once Po is known, one 
can combine equation (|23p with the deceleration law, e.g. 
ppn-2 _ cons t (equation I21|l . and solve by integration for 
the evolution of a given pulsar. This way one can produce 
isochrones, for example as a function of the braking index 
for pulsars with the same given value of Po. 
In Fig. [11] we have plotted pulsar isochrones for a pulsar 
with Po = 2.5ms and Po = 6.26 x 10 - 20 (e.g., correspond- 
ing to Bo = 1.15 x 10 s G for M = 1.4 Mq , sina = 1 and 
£o = 0) by varying the value of the braking index such that 
1 ^ n ^ 15. Our numerical calculations were confirmed by 
analytic calculations for the special cases where n — 2, 3 or 5 
(see small circles at the t = 12Gyr isochrone). The black 
dot indicates where the two 6 Gyr isochrones, calculated for 
n — 3 or Po = 6.26 x 10~ 20 , respectively, cross each other 
for a common solution and also in conjunction with the in- 
tersection of the evolutionary track calculated for the same 
values of n and Po. 



8.2 Characteristic versus true ages of MSPs 

Introducing the characteristic age of a pulsar, r = P/(2P) 
one finds the relation between r (the observable) and t: 

logr = log* + log - log (l - (^) n " 1 ) (24) 

for evolution with a constant braking index, n. The asymp- 
totic version of this relation (as t — ¥ oo and P 3> Po) is given 
by: 

{log t + log 2 for n — 5 
log t for n — 3 (25) 

logt-log2 forn = 2 

In Fig. [T2]we have plotted r as a function of t on log-scales 
for braking indices of n = 2, 3 and 5, respectively. Initially, 
when t is small and P ~ Po, r is always greater than t. After 
a certain timescale r can either remain larger or become 
smaller than t, depending on n. For n = 5 (n > 3) we 
always have r > t, for n — 3 we have r ^ t and only for 
n = 2 (n < 3) we have the possibility that r can be either 
greater or smaller than t. Having n = 2 corresponds to P 
being a constant, and solving for r = t yields P = 2Po and 
t — Po/ P- We notice that for a reasonable interval of braking 
index values 2 n ^ 5 the observed characteristic age will 
never deviate from the true age by more than a factor of two 
when P 3> Po- However, to reach P ^> Po may take several 
Hubble timescales if Bo is low. 

The characteristic age at birth (after recycling) is given 
by: tq — Po/(2Pq) and combining with equation (JS| we get: 
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where the constant k — 9.2 x 10 18 Gs -1 / 2 for sincto = 1 
and M = 1.4 M©. This expression verifies the well-known 
result that MSPs that are born relatively slowly spinning 
(large Po value) and/or have a relatively weak initial B- 
field (small Bo value) are also those MSPs which are born 
with characteristic ages which differ the most from their 
true ages as shown in Fig. [12] (see also Kizi ltan fc Thorsetj 
I2010T I. In other words, those pulsars with smallest values of 
Po for a given Bo or Po are those with to 2> t (which is 
hardly a surprise given that r = P/(2P)). The extent to 
which they continue evolving with r 2> t depends on their 
initial conditions (Po, Bo) and well as n. Notice, for recycled 
pulsars n is not measurable and for this reason its value 
remains unknown. 



8.3 Comparison with observations 

8.3.1 Kinematic corrections to P 

In a discussion of the ages and the evolution of millisec- 
ond pulsars in the PP-diagram kinematic corrections must 
be included when considering the observed values of P 
(Pobs). The kinematic corrections to the intrinsic P (Pint) 
are caus ed by acceleratio n due to proper motion of nearby 
pulsars (|Shklov skii 1970) and to vertical (az) and differen- 
tial rotational acceleration in our Galaxy (ogdr). The total 
corrections are given by: 



Pobs 



FU 
P 



+ 



Pshk 
~P~ 



az «GDR 

c c 



(27) 



and fo llowing lDamour fc Tavlorl l)l99ll ) and lWolszczan et al.1 
(2000) we can express these corrections as: 



-fobs 

~P~ 



az sin b 



(28) 



A. 

cRo 



cos I + (ri/Po)~cosZ 

l + (d/P ) 2 -2(d/P )cosZ 



where d is the distance to the pulsar, /i is the proper motion 
(related to the transverse velocity, v± = lid), az is the verti- 
cal component of Gala ctic acceleration (e.g. , from the Galac- 
tic potential model of iKuiiken fc Gilmorel (1 19891 )). I and b 
axe the Galactic coordinates of the pulsar, Po = 8.0 kpc and 
vo — 220kms _1 represent the distance to the Galactic cen- 
tre and the orbital velocity of the Sun, and c is the speed of 
light in vacuum. The corrections to P due to Galactic ver- 
tical and differential rotational accelerations are typically 
quite small (a few 10~ 22 <C Pint) and can be ignored, except 
in a few cases (see below). 



8.3.2 Evolutionary tracks and true age isochrones 

In order to investigate if we can understand the distribu- 
tion of MSPs in the PP-diagram we have traced the evo- 
lution of eight hypothetical, recycled MSPs with different 
birth locations. In each case we traced the evolution for 
2 ^ n ^ 5 and plotted isochrones similar to those intro- 
duced in Fig. 1111 The results are shown in Fig. 1 131 together 
with observed data. All the measured P values have been 
corrected for kinematic effects, as described above. If the 
transverse velocity of a given pulsar is unknown we used 
a value of 67kms _1 which we found to be the median 



value of the 49 measured velocities of binary pulsarfl In 
five cases (PSRs: J1231-1411, J1614-2230, J2229+2643, 
J1024-0719, J1801-1417, marked in squares in Fig.[T3J| we 
obtain P, n t < which is not physically possible. The reason 
is probably an overestimate of the pulsar distance. In those 
cases we have recalculated Pint assuming only half the dis- 
tance and included the corrections due to Galactic vertical 
and differential rotational acceleration. 
Three main conclusions can be drawn from this diagram: 

1) The overall distribution of observed pulsars follows nicely 
the banana-like shape of an isochrone with multiple choices 
for Po (or Po), see fat purple line. The chosen values of 
Po = 3.0 ms, n = 3 and t = 6 Gyr are just for illustrative 
purposes only and not an attempt for a best fit to the obser- 
vations. Fitting to one curve would not be a good idea given 
that MSPs are born with different initial spin periods, which 
depend on their accretion history, and given that the pulsars 
have different ages. The spread in the observed population 
is hinting that recycled pulsars are born at many different 
locations in the PP-diagram. 

2) The far majority of the recycled pulsars seem to have true 
ages between 3 and 12 Gyr, as expected since the population 
accumulates and the pulsars keep emitting radio waves for 
a Hubble time. 

3) Pulsars with small values of the period derivative P ~ 
10~ 21 hardly evolve at all in the diagram over a Hubble 
time. This is a trivial fact, but nevertheless important since 
it tells us that these pulsars were basically born with their 
currently observed values of P and P (first pointed out by 
ICamilo et al.lll994r i. In this respect, it is intere sting to no- 
tice P SR J1801— 3210 (recently discovered by iBates et al.1 
l201ll ) which must have been recycled with a relatively slow 
birth period, Po ~7 ms despite its low B-field < 10 8 G - see 
Fig. [5] for its location in the PP-diagram. 

One implication of the third conclusion listed above 
is that some radio MSPs must have been born with very 
small values of Bo — 1 x 10 7 G, if the inferred P-fields us- 
ing equation ((5]) are correct. Given their weak B-fields such 
sources would most likely not be able to channel the ac- 
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19981). One case 


radio pulsar J0034-0534 ( 


Bailes et al.lll994l 



transverse velocity of 146 km s 1 and correcting for the 
Shklovskii effect yields an intrinsic period derivative of only 
Pint = 5.36 x 10" 22 . These values result in B = 9 x 10 6 G 
(for M = 1.4 Mq and a = 90°) and a characteristic age, 
r = 55 Gyr. Once again we see that r is not a good mea- 
sure of the true age of a pulsar. If this pulsar has any 
gravitational wave emission (e 7^ 0) its derived surface B- 
field would be even smaller. We also notice a few pulsars 
which seem to have been recycled with M < 10 -3 A#Edd 
(approaching 10 _4 MEdd). Their progenitors could be the 
so-called weak LM XBs which are low-luminosity LMXBs 
(Ivan der Khslkooe). It is also quite possible that their po- 



4 If we add to this sample the measured velocities of 13 iso- 
lated pulsars which show strong signatures of being recycled 
(P < 100 ms and P < 10~ 16 ) the median velocity becomes 
69kms _1 , which would make the Shklovskii corrections slightly 
larger and thus Pint somewhat smaller. 
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Figure 13. Isochrones of eight hypothetical recycled pulsars born at the locations of the red stars. The isochrones were calculated 
for different values of the braking index, 2 ^ n ^ 5 (see Fig. 1111 for details). Also plotted are inferred B-field values (dashed lines) and 
characteristic ages, r (dotted lines). The thin grey lines are spin-up lines with M/Mrjdd = L 10 , 10 -2 , 10 — 3 and 10 — 4 (top to bottom, 
and assuming sina = <j> = u c = 1). In all calculations we assumed a pulsar mass of IAMq. It is noted that the majority of the observed 
population is found near the isochrones for t = 3 — 12Gyr, as expected. The fat, solid purple line indicates an example of a t = 6Gyr 
isochrone for pulsars with any value of Pq, but assuming Po = 3.0ms and n = 3. It is seen how the banana shape of such a type of an 
isochrone fits very well with the overall distribution of observed pulsars in the Galactic disk. Binary pulsars are marked with solid circles 
and isolated pulsars are marked with open circles. All values of P have been corrected for the Shklovskii effect using data from the ATNF 
Pulsar Catalogue. Pulsars in squares have been corrected for both a reduced distance and acceleration effects in the Galaxy. Pulsars with 
very small values of P are basically born on location in this diagram since their evolutionary timescale exceeds the Hubble time. Pulsars 
with relative small values of P and large values of P have the potential to constrain white dwarf cooling models - see text for discussion. 



sition in the PP- diagram was strongly shifted during the 
RLDP (see fig. 2 in lTaurisll2012T ) so that their average mass- 
accretion rate could have been larger prior to the RLDP. Al- 
ternatively, these MSPs demonstrate that <j> > 1 and ui c < 1 
which could move the spin-up line substantially downwards, 
as demonstrated in Fig. [S] 

It is obvious that a different sample of recycled pul- 
sar birth locations in Fig. [T3] would lead to somewhat 
different isochrones. However, we believe the qualitative 
interpretation is trustworthy and encourage further de- 
tailed statistical populati on analysis to be carried out (see 
iKiziltan fc Thorsettllioiol ). 

It is worth pointing out the difference in evolution- 
ary timescales for the progenitor systems of BMSPs with 
He WDs or CO WDs, respectively. The IMXBs, leading 
mainly to BMSPs with CO WD companions, evolve on nu- 
clear timescales of typically 100—300 Myr. Hence, BMSPs 
with CO WDs could have true ages from zero to up to about 
12 Gyr (assuming the first progenitor stars in our Galaxy 
formed ~ 1 Gyr after Big Bang) . The LMXBs, on the other 
hand, could easily have ages of 10 Gyr before they form, 



which is the typical timescale for a 1 Mq star to evolve into 
a sub-giant which fills its Roche-lobe in an orbit with P or b 
larger than a few days. One may then think that BMSPs 
with He WD companions must be much younger than BM- 
SPs with CO WD companions (the BMSP age being calcu- 
lated from after the RLDP). However, some LMXBs may 
have donor stars up to 2 Mq and these stars evolve on a 
timescale of only ~ 1.5 Gyr. Hence, BMSPs with He WD 
companions can have almost similar true (post recycling) 
ages as BMSPs with CO WD companions. 



8.3.3 PSR B1937+21 and PSR J0218+4232: two recently 
fully recycled pulsars? 

Although recycled pulsars with high values of P evolve 
fast across the PP-diagram, we should expect to ob- 
serve a few fully recycled pulsars with P ~ 10 -19 . The 
only two cases in the Galactic d isk known so far are 
PSR B1937+21 llBacker et al.lll982h and PSR 0218+4232 
|Navarro et al.1 1 19951) . see Fig. [5] The latter is a 2.32 ms 
pulsar (Pint = 7.66 x 10 -20 ) which not only has a small 
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characteristic age (r = 480 Myr) but also must have a 
fairly young true age. Even if we assume it evolved with a 
constant value of P (and disregard the possibility of n > 2 
which would have resulted in a larger initial value of Po 
and hence a younger age) we find t — 132 — 546 Myr when 
assuming Po is in the interval 1-2 ms. As expected, this 
system hosts a white dwarf companion star which is still 
hot enough to be detect ed. Interestingly en ough, all the 
cooling models applied bv lBassa et al.l |2003l ) to this white 
dwarf seem to indicate a somewhat larger age (see their 
fig. 4). This discrepancy has the intriguing consequence 
that either their applied cooling models are overestimating 
the age of the white dwarf, or this pulsar evolved with 
n < 2 and evolved upward in the PP-diagram. 

PSR B1937+21 is the first MSP discovered and is not 
only young, like PSR J0218+4232, but also isolated. This 
fast spinning pulsar (1.56 ms) has P n t = 1-05 x 10~ 19 
which yields a characteristic age of r = 235 Myr. We find 
upper limits (n = 2) for the true age to be in the inter- 
val t — 78 — 229 Myr when assuming Po is in the interval 
0.8-1.3 ms. MSPs (or any pulsar with a low P) are not 
believed to be formed directly from supernova explosions. 
They must interact with a companion star to spin up and 
decrease their B-field by a large amount. Thus it is interest- 
ing that this pulsar was able to evaporate its companion star 
within a timescale of (a few) 10 8 yr - unless it evolved from 
an ultra-compact X-ray binary possibly initially leaving a 
planet around it l|van Haaften et afll2012f ). 



8.3-4 PSR J1841+0130: a young mildly recycled pulsar 

PSR J 184 1+0 130 is the 29 ms binary pulsar discussed earlier 
in Section ^. 2. II It has an unusually high value of P, close to 
10~ 17 , which reveals that it is young and places it close to the 
spin-up line for M — MEdd (assuming sin a = <j> = ui c = 1), 
see Fig. The characteristic age of PSR J1841+0130 is 
58 Myr. An upper limit to its true age is 77 Myr, which 
is found by assuming n = 2 and Po = 10 ms (Po cannot 
have been much smaller than 10 ms under the above men- 
tioned assumptions regarding the spin- up line). If we assume 
n = 3, the upper limit for the true age is 51 Myr. The WD 
companion is therefore hot enough to yield spectroscopic in- 
formation that will reveal its true nature - i.e. whether this 
is a CO WD or a He WD, as discussed in Section 13.2.11 
According to the dispersion measure, DM — 125 cm -3 pc, 
of this pulsar its distance is about 3 kpc which might make 
such observations difficult. 



Table 3. Selected physical paramete rs of the binary millise cond 
pulsar J1614-2230 (data taken from iDemorest et al.ll2010h . 



Parameter 



value 



Pulsar mass 
White dwarf mass 
Orbital period 
Orbital eccentricity 
Pulsar spin period 
Period derivative 



1.97 ± 0.04 Mq 
0.500 ± 0.006 Mq 
8.6866194196(2) days 
1.30 ±0.04 x 10~ 6 
3.1508076534271 ms 
9.6216 x 10~ 21 



RLO. We find that AM cq = 0.11 Mq. This result is in fine 
accordance with our Case A calculation in Paper I where a 
total of 0.31 Mq is accreted by the neutron star. (Recall that 
AAicq is a minimum value for ideal, efficient spin up.) For 
the Case C scenario it is difficult to reconcile the required 
accretion of 0.11 Mq with the very small amount which is 
expected to be accreted by the neutron star evolving through 
a CE-phase (usually assumed to be <ti 10 -2 Mq). Further- 
more, in Paper I we argued that Case BB RLO (stable RLO 
from a naked helium star following a CE) was not an option 
for this system. 

The next check is to see if the spin-relaxation time 
scale was shorter than the mass-transfer timescale for 
PSR J1614— 2230. In order to estimate itorque we must have 
an idea of the B-field strength during the RLO. Since we 
have shown in Paper I that the mass-transfer rate in the fi- 
nal phase AB was near (slightly above) the Eddington limit, 
the value of B prior to phase AB must have been somewhat 
larger than its current estimated value of ~ 8.4 x 10 7 G, ac- 
cording to equation {5). Using B — 4 — 10 x 10 8 G we find 
ttorquc — 2 — 9 Myr. Since this timescale is shorter than the 
duration of RLO (~ 10 Myr for phase AB, see Paper I) we 
find that indeed it was possible for PSR J1614— 2230 to spin 
up to its equilibrium period. 

Based on both binary evolution considerations (Paper I) 
and the spin dynamics (this paper) we conclude that 
the Case A scenario is required to explain the exis- 
tence of PSR J1614-2230. From Fig. |9] we notice that 
PSR J1614— 2230 is an anomaly among the population of 
BMSPs with CO WD companions - it has an unusual rapid 
spin. This is explained by its formation via a stable and rel- 
atively long lasting IMXB Case A RLO (Paper I) which is, 
apparently, not a normal formation scenario for these sys- 
tems. 



9 SPINNING UP PSR J1614 2230 

Recent Shapiro dela y measurements of the radio millisecond 
pulsar J1614-2230 ijDemorest et al.ll2010l ) allowed a precise 
mass determination of this record high-mass neutron star 
and its white dwarf companion. A few key characteristic pa- 
rameters of the system are shown in Table [3] In Paper I 
we discussed the binary ev olution which le d to the forma- 
tion of this system (see also lLin et alj|201lh . For estimating 
the amount of necessary mass accreted in order to recycle 
PSR J1614— 2230 we can apply equation (|14[) and assuming, 
for example, an initial spin period of 2.0 ms following the 



9.1 Evolution of PSR J1614 2230 in the 
PP— diagram 

We have demonstrated in Paper I that PSR J1614-2230 
mainly accreted its mass during the final phase (AB) of 
mass transfer. As mentioned above, in this phase the mass- 
transfer rate was high enough that the accretion onto 
the neutron star was limited by the Eddington limit, i.e. 
M = Mfidd- A na t ural question to address (see also 
iBhalerao fc Kulkarnil (120111 )) is then how PSR J1614— 2230 
evolved to the present location in the PP-diagram which, 
at first sight, is even below the spin-up line expected for 
M = 10 -2 MEdd- However, one must bear in mind the de- 
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Figure 14. Simple models of how PSR J1614-2230 could have 
evolved from a birth near a spin-up line with M = M^dd to its 
currently observed location, shown by the red star, within 2 Gyr. 
The orange line indicates the location of the M = MEdd spin- 
up line if a = 10°, 4> = 1.4 and oj c = 0.25. PSR J1614-2230 is 
located close to this line which suggests it was born (recycled) 
close to it. The three black solid lines were calculated by using 
a = 90°, = 1 and ui c = 1. If PSR J1614-2230 was born near 
such a spin-up line with M = MEdd the need for subsequent 
spin-down after its birth would be severe and we have evolved a 
couple of toy models for this alternative: The green model repre- 
sents gravitational wave radiation and the blue model represents 
a phase of enhanced torque decay (see text). 



itational waves is caused by a time-varying quadrupole mo- 
ment which arises, for example, if the accretion onto the 
neutron star created no n-axisymmetric moments o f inertia 
relative to the spin axis <|Shapiro fc Teukolskvlll983T ) . In the 
second model we assumed an enhanced torque decay (n > 3) 
for a limited amount of time following the recycling phase. 
Such enhanced torque decay could be caused by either align- 
ment between the B-field axis and the spin axis of the pulsar 
(but not complete alignment, in which case no pulsar would 
be seen) or further decay of the surface B-field through some 
unspecified mechanism. In our calculation for the second 
model we assumed that effectively B or sin a decreased by 
a factor of about five on a decay timescale of 200 Myr. In 
Fig. [14] we show that the present position in PP-diagram 
can be reached from such toy models within 2 Gyr, which is 
in accordance with the cooling age of the CO WD inferred 
bv lBhalerao fc Kulkarnil (|201ll i. Notice, if PSR J1614-2230 
accreted 0.31 Mq, as suggested in Paper I, it would in prin- 
ciple be possible for it to spin up to a period of 0.92 ms 
(disregarding possible braking torque effects due to gravita- 
tional wave radiation. iBildstenl (1 19981 )) and thus the initial 
spin periods in our toy models are justified. However, we 
emphasize again that it is much less cumbersome to explain 
the current location of PSR J1614— 2230 in the PP-diagram 
simply by finetuning, in particular, <j> and lu c . Therefore, fu- 
ture observations of systems like PSR J1614— 2230 may help 
to constrain the disk-magnetosphere interactions. 



pendency of the parameters a, <j> and uj c when discussing the 
location of the spin- up line (cf. Section |4, 3 1 . 

If PSR J1614— 2230 has a magnetic inclination angle, 
a < 90°, and in case <j> « 1.4 and lu c rs 0.25, then the 
problem is solved since the M = MEdd spin-up line (see 
orange line in Fig. I14[) moves close to the current position 
of PSR J1614— 2230. Therefore we do not see any reason 
to question the recycling model of BMSPs, based on the 
ob servations of the PSR J161 4— 2230 system, as proposed 
bv lBhalerao k, Kulkarnil (|201ll ). However, if <j> « w c 1 and 
a fa 90° then we have to consider an alternative explanation 
which we shall now investigate. 

It is well-known that the rotational evolution of millisec- 
ond pulsars is rather poorly understood since the braking 
index only has been measured accurately for a few young, 
non-recycled pulsars. As briefly mentioned in Section[8]there 
is a variety of mechanisms which can influence the braking 
index of a pulsar. For example, the geometry of the B-field, 
the decay of the B-field, a changing magnetic inclination 
angle (i.e. between the spin axis and the B-field axis of 
the neutron star), the dynamics of the current flow in the 
pulsar magnetosphere and g ravi tational wave emission, e.g. 
iManchester fc Tavlorl l| 19771 ) and lContopoulos fc Spitkovskvl 
|2006t ). We have made a simple toy model for two of these 
mechanisms to show that PSR J1614— 2230 could in princi- 
ple have accreted with M — MEdd for sin a f» (j> f» oj c f» 1 
and subsequently evolve to its presently observed values of 
P and P (near the corresponding M = 10~ 2 MEdd spin-up 
line) after the radio pulsar turned on. In the first model we 
assumed pulsar spin-down as a result of gravitational wave 
radiation, which corresponds to n = 5. The emission of grav- 



10 NEUTRON STAR BIRTH MASSES 

In Paper I, we concluded that the neutron star in 
PSR J1614— 2230 must have been born in a supernova ex- 
plosion (SN) with a mass of 1.7± 0.15 Mq if it subsequently 
evolved through an X-ray phase with Case A RLO. Based 
on evolutionary considerations we also found an alternative 
scenario where the post-SN system could have evolved via 
Case C RLO leading to a CE. However, we also argued that 
the CE scenario was less probable given the difficulties in the 
required evolution from the ZAMS to the X-ray phase. In 
this paper we rule out the CE scenario completely in view of 
our analysis of the rotational dynamics. The in-spiral of the 
neutron star in a CE occurs on such a short timescale that 
no significant spin up is possible, and wind accretion from 
the proto WD cannot recycle the pulsar to a spin period of 
2—3 ms which requires accretion of 0.1 Mq. Recycling via 
a Case BB RLO, following CE evolution initiated while the 
progenitor of the donor star was on the RGB, is also not 
possible for PSR J1614— 2230 since its orbital period is so 
large that its post-CE orbit must have been wide too and 
thus there could not have been released enough orbital en- 
ergy during the in-spiral to eject the envelope. 
Does such a high birth mass of 1.7 Mq significantly exceed 
the neutron star birth masses in previously discovered ra- 
dio pulsar systems? Below we test if other pulsars may have 
been born in a SN with a similarly high mass. 

The interval of known radio pulsar masses ranges from 
1.17 M Q in the doub le neutron star binary PSR J1518+4909 
(3-<T upper limit, Ijanssen et alj 120081 ) to 1.97 Mq in 
PSR J1614— 2230, discussed in this paper. Prior to the mea- 
surement of PSR J1614— 2230 the highest mass inferred 
for the birth mass of a radio pulsar was only 1.44 Mq 
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Table 4. Measured masses of neutron stars in NS-WD systems in 
the Galactic disk, including PSR J1903+0327. The spin periods in 
the second column are in ms. The third, fourth and fifth column 
state the measured mass, the minimum accreted mass needed to 
spin up the pulsar using equation d 1 4 1) and the derived upper 
value for its birth mass, M^j 1 ^ = Mns — AAf oq , respectively. 
All masses are in quoted in Mq . All error bars are 1-<t values. 



Pulsar 


p . 

* spin 


M NS 


AM oq 


iyrbirth 


J0437-4715 1 


5.76 


1.76 ±0.20 


0.04 


1.72 


± 


0.20 


J0621+1002 2 


28.9 


1 70+ - 10 


0.005 


1.70 


± 


0.14 


J0751+1807 2 


3.48 


1.26 ±0.14 


0.07 


1.19 


± 


0.13 


J1012+5307 3 


5.26 


1.64 ±0.22 


0.05 


1.59 


± 


0.21 


J1141-6545 4 


394 


1.27 ±0.01 




1.27 


± 


0.01 


J1614-2230 5 


3.15 


1.97 ±0.04 


0.09 


1.88 


± 


0.04 


J1713+0747 6 


4.57 


-, co+O.OS 
1 " JO -0.06 


0.05 


1.48 


± 


0.07 


J1738+0333 7 


5.85 


1 46+ - 06 


0.04 


1.42 


± 


0.05 


J1802-2124 8 


12.6 


1.24 ±0.11 


0.01 


1.23 


± 


0.11 


B1855+09 9 


5.36 


1 -°'-0.11 


0.04 


1.53 


± 


0.11 


J1903+0327 10 


2.15 


1.667 ±0.021 


0.15 


1.52 


± 


0.02 


J1909-3744 11 


2.95 


1.44 ± 0.02 


0.09 


1.35 


± 


0.02 


B2303+46 12 


1066 


i no+0.06 
J-5O-0.10 




1.38 


± 


0.08 



References: (l)IVerbiest et ail 12003), (2) iNice et alj jjooj) , (3) 

Callanan et alj jl998T) , (4) lBhat etal] d2008h ■ ( 5) 

Demorest et alj 20ld) (6) ISplaver et alj 120051). ( 7) 

Antoniadis et al 120121), (8) iFerdman et alj 120101). (9) 

Nice et alj 120031), (lO)lFreire et alj l201ll) f 11) Ijacobv et alj 

120051) . (12) lThorsett fc Chakrabartvl ll999l) . 



(PSR 1913±16, see Paper I). Whereas most of the pulsars 
in double neutron star systems have their mass determined 
accurately, only a few of the ~ 160 binary pulsars with WD 
companions have measured masses. These are listed in Ta- 
bled At first sight, a handful of these are more massive than 
this previous limit of 1.44 Mq. However, in most of those 
cases the mass determinations have large uncertainties and 
in all cases also include mass accreted from the progenitor 
of their WD companion, i.e. these masses are pulsar masses 
after the recycling phase and not neutron star birth masses 
after the SN. 

To take this effect into account we have estimated upper 
limits for the birth mass of the neutron star in all systems 
by subtracting the minimum mass needed to spin up the 
pulsar using equation (|14[) . The results are shown in the 
fifth column of Table [4] Since we do not know the true age 
of the recycled pulsars, nor their braking index, we simply 
assumed in all cases that the presently observed pulsar spin 
period, P is related to their original equilibrium spin period, 
P eq prior to the RLDP via: P = V2P cq . For pulsars where 
the RLDP effect was small the initial spin period after recy- 
cling, Po — P cq and hence P ~ \/2Po, which corresponds to 
the case where the true age is half the characteristic age of 
a pulsar with braking index n = 3 (see Section [8] for further 
discussion on braking indices and true ages). 

It is possible that some BMSPs have P ~ Po (especially 
those BMSPs with small values of Bo). In this case the neu- 
tron star birth masses would be slightly larger than indicated 
here. However, recall that the limits in Table 2] are absolute 
upper limits obtained for idealized and efficient spin up. Fur- 
thermore, once an accreting pulsar reaches its equilibrium 



spin period it can still accrete material but the net effect may 
not be further spin up (depending on the mass-transfer rate, 
M). Therefore, it should be emphasized that the actual SN 
birth mass could in many cases be substantially lower than 
estimated in Table [4] For example, we have demonstrated in 
Paper I that PSR J1614-2230 could have accreted 0.31 Mq, 
which is substantially more than shown in Table [4] where the 
pulsar in this case was assumed to be born with a spin period 
of 3.15 ms/V2 = 2.23 ms, yielding AM eq = 0.09 Mq. 



10.1 PSR J0621 + 1002 

Considering the derived SN birth masses in Table [4] we 
notice at first sight in particular PSR J0437— 4715 and 
PSR J0621+1002 as further candidates for pulsars which 
were potentially born massive (~ 1.7 Mq), although the 
error bars are large. Given the uncertainty of the accre- 
tion efficiency (i.e. the possible difference between the ac- 
tual amount of mass accreted and AM eq ) we see that in 
particular PSR J0621±1002 could have been born massive 
with M^ma* ~ 1-7 M Q . PSR J0621+1002 (ICamilo et all 



Il99rj ) has an orbital period of 8.3 days and a CO WD com- 
panion of mass 0.67 Mq which is somewhat equivalent to 
PSR J1614— 2230. However, its relatively slow spin period 
of 28.9 ms and large eccentricity of 2.5 x 10 -3 indicate that 
PSR J0621±1002 did not evolve through an X-ray phase 
with stable, and relatively long Case A RLO (unlike the 
situation for PSR J1614— 2230). It seems much more evi- 
dent that PSR J0621+1002 evolved through a rapid, less 
efficient, mass-tran sfer phase - such as early Case B RLO 
l|Tauris et al.ll20o"oh . In this case the actual amount of accre- 
tion must have been at most a few 10 -2 Mq. Therefore, its 
measured mass is very close to its birth mass and if future 
mass measurements of PSR J0621+1002 yield smaller error 
bars, still centered on its current mass estimate, we conclude 
that this pulsar belongs to the same class of neutron stars as 
PSR J1614-2230, Vela X-l and possibly the Black-Widow 
pulsar which were all born massive (see Paper I). However, 
the error bar is quite large and PSR J0621+1002 could turn 
out to have an ordinary mass of about 1.4 Mq (the 2-a lower 
mass limit). 



10.1.1 PSR J0437-4715 

PSR J0437— 4715 is also a potential candidate for a pul- 
sar which could have been born as a neutron star with a 
mass of ~ 1.7 Mq. However, this pulsar has a fairly rapid 
spin period (like a few other relatively high-mass neutron 
star candidates in Table [4} which indicates a long spin-up 
phase. Hence, it may have accreted significantly more than 
suggested by AM oq . In fact, as mentioned before, in the 
case of PSR J1614— 2230 the pulsar might have accreted 
three times more than suggested by AM oq , see fig. 7 in 
iTauris et all < |201lT >. The problem is that as soon as the pul- 
sar spin reaches P eq (which requires accumulation of mass 
AM oq ) it may remain in near-spin equilibrium and keep ac- 
creting on and off while spinning up / down depending on the 
sign of the accretion torque. The situation is quite different 
for PSR J0621+1002 which formed through a short lived X- 
ray phase excluding evolution at near-spin equilibrium for 
a substantial amount of time. We therefore conclude that 
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PSR J0437-4715 could originally have been born in a SN 
with a typical mass of ~ 1.4 Mg, even if its present mass is 
indeed ~1.76M Q . 



11 SUMMARY 

We have investigated in detail the recycling process of pul- 
sars with respect to accretion, spin up and the Roche-lobe 
decoupling phase, and discussed the implications for their 
spin periods, masses and true ages. In particular, we have 
discussed the concept of a spin-up line in the PP-diagram 
and emphasize that such a line cannot be uniquely defined. 
Besides from the poorly known disk-magnetosphere physics, 
which introduces large uncertainties, for each individual pul- 
sar the equilibrium spin period, P cq also depends on its mag- 
netic inclination angle, a as well as its accretion history (M) 
and its B-field st rength. Furthermore, we have applied the 
ISpitkovskvl HH) spin-down torque on radio pulsars which 
significantly changes the location of the spin-up lines com- 
pared to using the vacuum magnetic dipole model, especially 
for small values of a. 

We have derived a simple analytical expression (equa- 
tion I14[) to evaluate the amount of mass needed to be ac- 
creted to spin up a pulsar to any given equilibrium s p in pe - 
riod, P cq . Our result resembles that of lAlpar et al.l (|l982h 
and approximately yields the same v alues (within a fac- 
tor of two) as the expression derived bv lLipunov &: Postnovl 
l| 19841 ). Using our formula we find, for example, that BM- 
SPs with He WDs and P eq — 2 ms must have accreted at 
least 0.10 M Q , whereas typical BSMPs with CO WDs and 
P C q ~ 20 ms only needed to accrete 0.005 Mq. 

Applying equation (|14fl enables us to explain the dif- 
ference in spin distributions between BMSPs with He and 
CO WDs, respectively. The BMSPs with He WDs often 
evolved via an X-ray phase with stable RLO on a long 
timescale, allowing sufficient material to be accreted by the 
neutron star to spin it up efficiently to a short period, 
whereas the BMSPs with CO WDs most often have slow 
spins and evolved from IMXBs which had a short phase of 
mass transfer via early Case B RLO or Case C RLO - the 
latter leading to a CE-evolution often followed by Case BB 
RLO. The only exception known so far is PSR J1614-2230 
since this system produced a CO WD orbiting a fully recy- 
cled MSP and thus it must have evolved via Case A RLO of 
an IMXB. 

It is not possible to recycle pulsars to become MSPs 
via wind accretion from 1.1 — 2.2 M© post-CE helium stars. 
However, we have demonstrated that Case BB RLO from 
such helium stars can spin up pulsars to at least ~ 11 ms. 
Further studies of these systems are needed. 

There is an increasing number of recycled pulsars with 
WD companions which seem to fall outside the two main 
populations of BMSPs with He and CO WDs, respectively. 
These peculiar systems possibly have He WDs and always 
exhibit slow spin periods between 10 and 100 ms. We suggest 
that these systems with P orD . ^ 1 day may have formed via 
Case A RLO of IMXBs. We plan further studies on these 
binaries. 

The Roche-lobe decoupling ph ase (RLDP), at the ter- 
minal stage of the mass transfer ijTaurid [2012 ). has been 
analysed and we have shown that while the RLDP effect 



is important in LMXBs - leading to significant loss of rota- 
tional energy of the recycled pulsars as well as characteristic 
ages at birth which may exceed the age of the Universe - 
it is not significant in IMXB systems where the duration of 
the RLDP is short. 

In order to track the evolution of pulsars in the PP- 
diagram we have introduced two types of true age isochrones 
- one which matches well with the banana shape of the ob- 
served distribution of known MSPs. We encourage further 
statistical population studies to better understand the for- 
mation a nd evolution of radio MS Ps in the PP-diagram 
(see e.g. iKiziltan fc Thorsettl l2010h . The discrepancy be- 
tween true ages and characteristic spin-down ages of recycled 
pulsars has been discussed and we confirm that the latter 
values are completely untrustworthy as true age indicators, 
leaving WD cooling ages as the only valid, although not ac- 
curate, measuring scale (|Taurisll2012T ). 

In the combined study presented here and in Paper I we 
have investigated the recycling of PSR J1614— 2230 by de- 
tailed modelling of the mass exchanging X-ray phase of the 
progenitor system. Given the rapid spin of PSR J1614— 2230 
(3.15 ms) we argue that it is highly unlikely that it evolved 
through a CE, leaving Case A RLO in an IMXB as the 
only viable formation channel. We confirm the conclusion 
from Paper I that the neutron star in PSR J1614— 2230 was 
born massive (1.70 ± 0.15 Mq). We have demonstrated that 
PSR J1614-2230 could have been spun-up at M = M E d d 
and subsequently evolve to its current position in the PP- 
diagram within 2 Gyr (the estimated cooling age of its white 
dwarf companion) . 

Besides PSR J1614-2230, Vela X-l and possibly 
the Black-Widow pulsar, we have argued that also 
PSR J0621+1002 could belong to the same class of neu- 
tron stars born massive (> 1.7 Mq). The formation of such 
massive neutron stars in supernovae is in agreement with 
some supernova ex plosion models (e.g. IZhang et al.l 120081 ; 
lUgliano et aill2012l ). 
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APPENDIX A: IDENTIFYING THE NATURE 
OF A PULSAR COMPANION STAR 

For stellar evolutionary purposes, and for observers to judge 
the possibilities to detect a companion, it would be useful to 
have a tool to predict the most likely nature of the compan- 
ion star of a given binary pulsar. As discussed in Section [3] 
radio pulsars are presently found with the following com- 
panions: 

• Main sequence stars - MS 

• Neutron stars - NS 

• CO/ONeMg white dwarfs - CO 

• He white dwarfs - He 

• Ultra-light companions (or planets) - UL 

Although still not detected, radio pulsars are also expected 
to exist in systems with either: 

• Helium stars - HS 

• (sub) Giant stars - GS 

• Black Holes - BE 

As explained in Section [3] pulsar systems found in globular 
clusters are not useful as probes of stellar evolution due to 
perturbations in their dense environment. We shall denote 
the companions of such pulsars by: 

• Globular cluster companions - GC 

The usual characteristic parameters measured of 
an observed binary pulsar are: the spin period (P), the 
period derivative (P), the orbital period (P or b) and the 
projected semi- major axis of the pulsar (01). From the 
latter two parameters one can calculate the mass func- 
tion (/) which is needed, but not sufficient, to estimate 
the mass of the companion star (M2). Furthermore, in 
many cases it is also possible to estimate the eccentric- 
ity (ecc). Based on these parameters we state simple 
conditions in Table I All which can be used to identify 
the most likely nature of pulsar companions at the 
95 per cent level of confidence (based on statistics from 
current available data). These conditions can be directly 
applied to the online, public available ATNF Pulsar Cat- 
alogue, http://www.at nf.csiro.au/research / pulsar /psrcat / 
((Manchester ct al. 2005) in order to select a subsample of 
pulsars with certain companion stars. 

One example of the few pulsars where the identification 
fails is PSR J1614— 2230 given its unusual combination of 
being fully recycled and having a CO WD companion, as 
discussed in this paper. The identification conditions can in 
principle be optimised further to increase the level of confi- 
dence even more. However, this would not only require these 
conditions to be updated frequently with new discoveries of 
somewhat peculiar binary pulsars but also lead to unde- 
sirable complexity in the conditions. The main uncertainty 
factor is the (most often) unknown orbital inclination angle 
of the binary system. The transition between systems with 
UL and He companions and, in particular, between He and 
CO companions is difficult to define satisfactory. Further- 
more, WD+NS systems where the last formed compact ob- 
ject is the neutron star are very difficult to distinguish from 
NS+NS systems, i/based solely on data from the last formed 
(non-recycled) NS. Finally, we note that the long-sought af- 
ter radio pulsar in a binary with a black hole companion 



(NS+BH or BH+NS) will most likely be mildly recycled 
(i.e. NS+BH) due to its much longer lifetime as an observ- 
able pulsar after accretion. In this case it may be difficult 
to distinguish it from a system where the pulsar has a NS 
companion orbiting with a small inclination angle. On the 
other hand, if the pulsar is not recycled then such a BH+NS 
system could resemble a pulsar with an (undetected) MS 
companion. 
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Table Al. The most likely nature of a radio pul sar companion star base d on observable characteristics. The selection conditions can be 
directly applied to the ATNF Pulsar Catalogue (Manchester ct al. 2 0051 ). See also Section [3] for further discussions. 



Companion type Conditions 



MS 


M* 


> 0.5 M Q 


and 


P orb > 50 d 








NS 








P orb < 50 d 


and 


ecc > 0.05 




CO 


M* 


> 0.335 M 


and 


^orb < 75 d 


and 


ecc < 0.05 and 


P > 8 ms 


He 


M* 


> 0.08 M Q 


and 


{(Porb < 75 d 


and 


Ml < 0.335 M Q ) 


or (P orb > 75 d and M% < 0.46 M )} 


UL** 


M* 


< 0.08 M 













* The median companion mass M2 is calculated for an orbital inclination angle i = 60° and an assumed pulsar mass Mns = 1.35 Mq. 
** Many pulsars with unmeasured values of a\ are also expected to host an ultra-light companion if P or b < 2 days and P < 8 ms. 
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